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ABSTRACT 

An  experimental  investigation  was  conducted  to  study  convective  heat  transfer  in 
straight  and  curved  rectangular  channels  of  high  aspect  ratio  that  approximate  plates 
of  infinite  extent.  Experiments  were  performed  at  steady  state  in  the  turbulent  flow 
regime  with  one  wall  held  at  a  constant  heat  fiux  and  the  opposite  wall  essentially 
adiabatic.  The  effect  of  curvature  induced  secondary  flow  on  heat  transfer  on  the 
concave  and  convex  walls  was  observed  by  comparing  Nusselt  numbers  for  four 
different  configurations  at  several  different  Reynolds  numbers.  Significant  heat  transfer 
enhancement  was  observed  on  the  concave  wall.  Correlations  for  Xusselt  number  as  a 
function  of  Reynolds  number  were  calculated  for  the  cases  studied. 
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I.  INTRODUCTION 

A.       BACKGROUND 

Lord  Rayleigh  [Ref.  1]  was  seemingly  the  first  to  investigate  the  mechanics  of 
Hows  with  curved  streamlines.  His  analysis  showed  that  inviscid  curved  streamline 
flows  are  stable  providing  that  circulation  increases  with  radial  distance  (for  flow 
between  concentric  cylinders  this  condition  is  met  when  the  outer  cylinder  spins  more 
quickly  than  the  inner  cylinder). 

Later.  Taylor  [Ref  2]  studied  viscous  flow  between  concentric  rotating  cylinders 
and  showed  that  when  the  inner  cylinder  is  rotated  while  holding  the  outer  cyUnder 
stationan.'.  the  onset  of  instability  is  predicted  by  the  dimensionless  parameter: 

Ud        /~ 
Ta  (Tavlor  number)  =    —      v —        >41.3,  (eqn  1.1) 

V  r. 

where  U  is  the  tangential  velocity  of  the  inner  cylinder,  d  is  the  gap  between  cylinders, 
r.  is  the  radius  o[  curvature  of  the  inner  cylinder,  and  v  is  kinematic  viscosity.  The 
instability  is  manifest  in  the  formation  of  regularly  spaced,  counter-rotating  vortices 
whose  axes  of  rotation  lie  in  the  circumferential  (main  flow)  direction  (see  Figure  1.1). 
Furthermore,  Taylor  concluded  that  these  vortices  form  in  the  laminar  flow  regime  but 
do  not  induce  turbulent  main  flow  [Ref  3]. 

The  flow  of  a  viscous  fluid  in  a  curved  channel  was  studied  by  Dean  [Ref  4].  He 
developed  an  analytical  solution  for  flow  in  a  curved  channel  approximated  by 
concentric  cylinders  where  the  radius  of  the  inner  cylinder  was  much  greater  than  the 
gap  between  cylinders.  Instead  of  rotating  the  inner  cylinder  Dean  imposed  a  pressure 
gradient  along  the  channel  wall  and  predicted  that  flow  instability  would  occur  when  a 
critical  value  o[  a  dimensionless  parameter  similiar  to  the  Taylor  number  was  exceeded: 


Ud        fd~ 

De  (Dean  number)  =    —    v —        >  36.  (eqn  1.2) 

V  r. 

1 


where  d  represents  the  channel  half-width. 

12 


Figure  1.1     Taylor  Vortices  Between  Rotating  Cylinders. 
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Gortler  [Ref.  5]  went  on  to  investigate  viscous  flow  in  a  curved  channel  and 
concluded  that  the  flow  instability  predicted  by  Dean  would  take  a  form  analagous  to 
the  Taylor  vortices  of  concentric  cylinder  flow  (see  Figure  1.2).  Regularly  spaced 
counter  rotating  vortices  form  along  axes  of  rotation  parallel  to  the  main  flow 
direction.  Yih  and  Sangster  [Ref.  6]  concluded  that  the  centrifugal  force  coupled  with 
the  velocity  gradient  near  the  concave  wall  caused  instability  due  to  a  stratification  in 
specific  weight  that  led  to  vortex  formation. 

The  secondan.-  flows  predicted  by  Taylor  and  Gortler  (Taylor-Gortler  vortices) 
have  been  confirmed  experimentally  using  various  flow  visualization  techniques,  and 
hot  wire  anemometr}'  [Ref.  7].  The  presence  of  these  vortices  in  the  laniinar  flow 
regime  led  to  speculation  ot^  the  eflect  of  this  secondan.'  flow  on  heat  transfer  over 
curved  surfaces.  Extensive  research  has  been  done  which  yields  good  experimental 
agreement  with  analytical  solutions  for  laminar  flow  in  curved  channels  of  rectangular 
cross  section  [Refs.  8.9,10],  and  for  turbulent  flow  in  a  curved  rectangular  channel  of 
low  aspect  ratio  (ratio  of  the  spanwise  dimension  to  the  channel  height) 
[Refs.  11,12,13].  Research  in  these  areas  indicates  that  there  is  enhancement  of  heat 
transfer  in  these  flow  regimes  when  Taylor-Gortler  vortices  form  near  the  concave  wall. 
However,  information  relating  to  the  fluid  mechanics  and  heat  transfer  in  the  transition 
and  turbulent  flow  regimes  for  channels  of  high  aspect  ratio  is  ver>'  limited. 

B.  APPLICATIONS 

Studying  heat  transfer  over  a  curved  surface  is  just  one  aspect  of  the  complex 
fluid  dynamics  and  heat  transfer  problems  that  exist  in  an  operating  gas  turbine.  The 
high  temperatures  characteristic  of  a  combustion  engine  constrain  gas  turbine 
performance  because  of  metalurgical  limitations  of  turbine  blades  which  in  turn 
restricts  the  efllciency  of  the  turbine.  Information  gathered  in  this  type  of 
experimentation  could  lead  to  a  better  understanding  of  the  problem  and  facilitate  the 
design  of  gas  turbines  of  higher  efficiency  and  power  to  weight  ratio. 

C.  INTENT  OF  THIS  STUDY 

The  purpose  of  this  investigation  was  to  study  the  eflect  of  curvature  on  heat 
transfer  in  a  curved  rectangular  channel  ol'  large  aspect  ratio  by  gathering  extensive 
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Figure  1.2    Taylor-Gortler  Vortices  in  a  Curved  Channel, 
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data  relating  Nusselt  number  to  Reynolds  number  for  flow  over  walls  of  constant  heat 
flux  of  straight,  convex  and  concave  geometry.  This  information  will  provide  a  firm 
data  base  for  validating  computer  models  of  the  complex  flow  and  heat  transfer  in 
curved  rectangular  channels  similar  in  geometry. 

Data  for  hydraulically  fully  developed  flow  in  a  straight  rectangular  channel 
heated  either  from  above  or  below  was  taken  to  provide  a  basis  for  comparison  with 
the  flow  over  a  heated  convex  or  concave  wall.  Empirical  correlations  similiar  to  the 
familiar  Dittus-Boelter  correlation  for  Nusselt  number  as  a  function  of  Reynolds 
number  were  calculated  from  the  data  so  that  the  heat  transfer  enhancement  could  be 
described  quantitatively  and  compared  to  previous  data. 

A  secondan.'  goal  was  to  determine  the  relative  eflect  of  buoyancy  on  heat 
transfer  in  a  curved  channel.  The  eflect  of  buoyancy  on  heat  transfer  in  forced  flow 
must  be  considered  if  the  eflects  of  free  and  forced  convection  are  comparable  because 
the  combined  heat  transfer  is  then  dependent  on  the  direction  of  free  convection  in 
comparison  to  that  of  forced  convection. 


II.  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

A  detailed  description  of  the  materials  and  construction  methods  used  to  build 
the  experimental  channel  and  associated  instrumentation  is  contained  in  G.  Galyo's 
thesis  [Ref.  14:  pp.  25-43].  Some  minor  changes  had  to  be  made  to  the  existing 
hardware  and  computer  software  and  are  described  herein. 

A.       DESCRIPTION  OF  THE  APPARATUS 

The  channel  was  built  of  plexiglass,  having  a  channel  width  of  25.4  cm  and  a 
height  of  0.635  cm,  yielding  an  aspect  ratio  of  40:1.  The  channel  was  rigidly  supported 
on  a  wooden  frame  and  insulated  over  its  entire  length  with  three  layers  of  half  inch 
foam  insulation.  The  curved  section  subtends  180°  of  arc,  and  the  convex  wall  (inner 
curved  wall)  had  a  radius  of  curvature  of  29.7  cm  (see  Figure  2.1).  Air  was  pulled 
through  the  channel  by  a  compressor  joined  to  the  duct  by  two  inch  PVC  pipe  through 
a  throttle  valve. 

Four  30.48  cm  x  25.4  cm  wire  wound  rubber  heaters  laminated  to  copper  plates 
were  mounted  on  the  inside  walls  of  the  channel  to  accomplish  heating  of  the  air.  Two 
heaters  were  mounted  in  the  straight  test  section,  one  on  the  upper  wall  and  one  on  the 
lower,  at  a  sufficient  distance  from  the  channel  entrance  to  insure  that  flow  was 
hydrodynamically  fully  developed  prior  to  flowing  over  the  straight  heated  sections. 
The  two  remaining  heaters  were  installed  in  the  curved  test  section,  one  on  the  concave 
wall  and  one  on  the  convex  wall,  with  the  leading  edges  of  the  heaters  beginning  at  an 
angle  of  101.5"  from  the  curved  section  entrance;  this  allowed  vortices  to  form  prior  to 
flowing  over  the  curved  heated  sections.  The  heaters  were  constant  power  density  type 
capable  of  600  watts  each,  and  could  be  energized  individually  or  in  any  combination 
of  pairs  by  use  of  two  power  supplies. 

Thirty  copper-constantan  thermocouples  were  imbedded  in  each  of  the  copper 
plates  at  the  copper  plate-rubber  heater  interface  so  that  local  temperatures  could  be 
measured  (see  Figure  2.2).  Four  thermocouples  were  mounted  laterally  across  the 
channel  entrance  to  measure  inlet  temperature  (see  Figure  2.1).  Similiarly.  four 
thermocouples  had  been  mounted  at  the  exit  of  the  channel;  however,  experimentation 
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Figure  2.1     Experimental  Rectangular  Channel. 
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Figure  2.2    Thermocouple  Placement  in  the  Heated  Plates. 
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showed  that  there  was  insufiicient  mixing  of  the  air  up  to  the  exit  which  resulted  in 
erratic  outlet  temperature  measurements  [Ref.  14:  p.  67].  Instead,  eighth  inch 
honeycombing  was  installed  at  the  exit  of  the  channel  to  promote  mixing,  and  the  four 
exit  thermocouples  were  placed  in  the  PVC  pipe  near  the  end  of  the  channel  outlet 
transition  bell  (see  Figure  2.1).  Two  layers  of  half  inch  foam  insulation  were  applied  to 
the  PVC  to  a  distance  o[  approximately  fourty-eight  centimeters  past  these 
thermocouples  to  minimize  heat  loss  to  the  surroundings  prior  to  measuring  outlet 
temperature.  One  thermocouple  was  placed  downstream  of  the  flow  measuring  orifice 
and  used  to  calculate  the  density  of  the  air  flowing  through  the  orifice.  Three  groups 
of  five  thermocouples  each  were  connected  in  parallel  between  the  layers  of  insulation 
surrounding  the  heated  sections  so  that  losses  to  the  surroundings  through  the 
insulation  could  be  estimated. 

Any  combination  of  eighty  voltages  from  the  thermocouples  and  power  supplies 
could  be  read  by  the  Hewlett  Packard  3054A  automatic  data  acquisition  system  used 
for  this  research.  Coupled  to  the  data  acquisition  system  were  a  Hewlett  Packard  9826 
computer  and  a  printer  which  provided  control  of  the  data  acquisition  system  and 
facilitated  recording,  reducing,  and  printing  of  data. 

A  flow  measuring  orifice  was  installed  in  accordance  with  the  AS  ME  Power  Test 
Code  [Ref  15],  and  water  and  mercurv'  manometers  were  used  to  measure  the  pressure 
upstream  of  the  orifice  and  dilTerential  pressure  across  the  orifice  so  that  mass  flowrate 
could  be  calculated. 

B.       EXPERIMENTAL  PROCEDURE 

The  procedures  followed  for  testing  were  the  same  irrespective  of  the  specific  test 
configuration.    The  procedures  are: 

1.  Zero  the  manometers,  shut  the  system  throttle  valve,  then  start  the  compressor. 

2.  Use  the  throttle  valve  to  adjust  the  diflerential  pressure  across  the  How 
measuring  orifice  to  that  corresponding  to  the  desired  flowrate  or  Reynolds 
number.  Determining  the  approximate  diflerential  pressure  corresponding  to  a 
specific  Reynolds  number  takes  some  experience  and  the  use  of  the 
approximation: 
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Re 

(AP)2  --  [— V{AP)i  (eqn2.1) 

Re, 

3.  Energize  the  heater  and  adjust  the  voltage  for  the  approximate  value  of 
(T  ,  -  T.  )  desired.  In  these  experiments  a  value  of  six  to  twelve  degrees 
centigrade  was  used  for  the  temperature  difierential,  the  objective  being  to 
maximize  (T  -T.  )  while  limiting  the  highest  local  plate  temperature  to  a 
maximum  of  about  lOO'C.  Maximizing  this  temperature  differential  lessens  the 
uncertainty  of  the  calculated  Nusselt  number. 

4.  Experience  showed  that  the  channel  reached  steady  state  about  four  hours  after 
first  energizing  the  heater,  and  approximately  one  hour  to  steady  out  between 
consecutive  data  runs  at  different  flowrates.  When  all  of  the  following  criteria 
were  met.  steady  state  was  considered  to  have  been  achieved: 

a.  Change  in  the  quantitv  (T     -T.  )  of  less  than  two  percent   over  a  ten 
minut^  period.  '        ^"^     '" 

b.  Chanee  in  the  heated  wall  temperature  of  less  than  one  percent  over  a  ten 
minut'e  period. 

c.  Chanse  in  heated  wall  insulation  temperatures  of  less  than  one  percent  over 
a  ten  niinute  period. 

5.  When  the  above  criteria  for  steady  state  are  met,  run  the  data  gathering  and 
reducing  program  (see  Appendix  A)  several  times  in  succession  to  verify  that 
steady  state  has  been  reached.  Three  to  six  consecutive  data  collection  runs  at 
steady  state  meeting  the  criteria  of  (4)  above  are  required  to  determine  an 
average  Nusselt  number  versus  Reynolds  number  relationship.  Generally  the 
run-to-run  variation  in  Nusselt  number  and  Reynolds  number  were  much  less  at 
lower  fiowrates. 

6.  When  a  sufficient  number  oC  consistent  data  runs  have  been  recorded,  throttle 
to  adjust  flowrate  to  the  next  desired  Reynolds  number  and  adjust  the  heater 
voltage  to  achieve  the  desired  plate  temperature. 

The  atmospheric  pressure,  orifice  dilferential  pressure  and  pressure  upstream  o[ 
the  orifice  had  to  be  read  from  a  barometer  and  two  manometers  and  entered  into  the 
data  reduction  program  interactively;  these  parameters  were  not  read  by  the  automatic 
data  acquisition  system. 
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During  this  study  it  was  noted  that  small  (one  degree)  ambient  temperature 
perturbations  had  a  noticeable  etlect  on  temperatures  measured  throughout  the  test 
system;  therefore,  close  control  of  ambient  temperature  oscillations  was  required  so 
that  temperature  changes  in  the  test  system  could  be  minimized  and  consistent  data 
measured  from  run  to  run. 

The  initial  configuration  of  the  test  channel  was  such  that  buoyancy  induced  flow 
opposed  forced  flow  in  the  curved  test  section.  After  data  was  taken  over  the  full 
range  of  fiow  that  the  compressor  was  capable  of.  the  rig  was  inverted  and  data  was 
taken  with  buoyancy  induced  flow  aiding  forced  fiow.  The  purpose  of  this  aspect  of 
the  study  was  to  determine  whether  natural  convection  contributed  significantly  to  the 
total  heat  transfer  in  the  curved  test  sections. 


III.  DATA  REDUCTION 

Each  heater  plate  provided  a  constant  surface  heat  flux.  The  heater  plates  were 
energized  by  regulated  D.C.  power  supplies  through  stable  precision  resistors  of  low 
resistance  value  (see  Figure  3.1).  The  voltage  drop  across  both  the  heater  element  and 
the  precision  resistor  were  read  and  recorded  by  the  data  acquisition  system.  Then  the 
total  power  supplied  by  the  heater  was  determined  from  the  relation: 

V.    X  V 

Qp  =    ^    ^    P*"  (eqn3.I) 

where  R     was  a  known  value  of  resistance  included  as  a  constant  in  the  data  reduction 
pr 

program.  The  insulated,  unheated  walls  were  treated  as  adiabatic  because  the  heat  loss 
through  them  was  negligible  compared  to  the  total  heat  input  from  the  heater  plate. 
Therefore,  the  bulk  temperature  varied  linearly  from  the  leading  to  the  trailing  edge  of 
the  heated  plate. 

Thermocouple  voltages  were  converted  to  centigrade  temperature  using  individual 
thermocouple  calibration  curves.  All  temperatures  were  recorded  then  combined  to 
calculate  specific  average  temperatures  desired. 

The  actual  heat  loss  to  the  surroundings  was  difficult  to  account  for.  Typically 
the  estimated  total  heat  loss  by  conduction  through  the  insulation  and  by  radiation 
was  only  about  ten  percent  of  the  actual  total  heat  loss.  The  major  difficulty  in 
making  an  accurate  prediction  of  the  heat  loss  was  due  to  not  being  able  to  determine 
the  area  of  the  channel  over  which  heat  was  lost  to  the  surroundings.  For  the  purpose 
of  estimating  the  heat  loss,  the  area  of  the  heated  plate  was  used  as  the  effective  area 
of  loss;  however,  the  plexiglass  section  of  the  channel  adjacent  to  the  heater  plates  acts 
as  extended  heat  transfer  surfaces,  and  the  other  copper  heater  surfaces  in  the  channel, 
which  have  a  low  resistance  to  thermal  conduction,  act  like  heat  sinks.  Therefore, 
instead  of  assuming  that  the  energy  transferred  to  the  air  was  equal  to  the  energy  input 
from  the  plate  minus  the  estimated  energy  lost  to  the  surroundings  over  the  area  of  the 
plate,  the  energy  actually  transferred  to  the  working  fluid  was  calculated  as  follows: 


23 


vWaa 


Rubber   Heater 


Precision  Resistor 


Power  Supply 


Figure  3.1     Po\^er  Measuring  Circuit. 


24 


Qair  "  mCpAT  (eqn  3.2) 

where  AT  =  T  ,    -  T. 

out        in 

The  corresponding  plate  average  heat  transfer  coefTicient  and  Nusselt  number  based  on 
the  hydraulic  diameter  of  the  duct  were  then  calculated  as  follows: 


Q 

h  =  ^^^ (eqn  3.3) 

—  hD.  . 

^%d  =  -^  (eqn  3.4) 

air 

where  D^^,  the  hydraulic  diameter  of  the  channel,  is  a  geometric  constant, 

4  X  Cross  Sectional  Area  of  the  Duct 

D.  .  =  (eqn  3.5) 

"'^  Wetted  Permieter  of  the  Duct 


Data  on  atmospheric  pressure,  orifice  differential  pressure,  and  orifice  upstream 
pressure  was  used  with  the  flow  measuring  orifice  characteristics  and  air  density  based 
on  air  temperature  at  the  orifice  to  calculate  mass  flowrate,  m.  Reynolds  number 
based  on  hydraulic  diameter  was  then  found  by  the  relation: 

mD.  . 
c^air 

where  ji^j^  was  the  viscosity  of  the  working  fluid  at  the  average  bulk  temperature 
measured  in  the  test  channel.  Siniiliarly.  Reynolds  numbers  based  on  the  pipe  diameter 
and  channel  height  were  also  determined.  For  the  concave  and  convex  sections  only, 
the  dimensionless  parameter.  Dean  number,  was  calculated  from: 


m  D  /2         /D  /2 

De  =  2 —      V-^ —  (eqn  3.7) 

A  u  .  r. 


25 


IV.  EXPERIMENTAL  RESULTS 

Several  sets  of  data  were  taken  over  the  approximate  range  of  Reynolds  numbers 
from  6000  to  20,000.  A  computer  program  (see  Appendix  A)  was  used  to  reduce  the 
raw  data  and  calculate  the  quantities  of  interest,  particularly  plate  average  Nusselt 
numbers  for  corresponding  Reynolds  numbers,  both  based  on  the  hydraulic  diameter  of 
the  channel.  Three  to  six  separate  data  sets  at  each  Reynolds  number  were  taken,  then 
averaged  to  give  the  value  of  the  quantities  of  interest. 

After  extensive  data  had  been  taken  for  the  channel  configured  such  that  free 
convection  opposed  the  direction  of  forced  convection,  the  channel  was  inverted  and 
several  sets  of  data  were  taken  for  the  convex  and  concave  plates  with  the  directions  of 
free  and  forced  convection  coincident  so  that  the  etTect  of  buoyant  flow  could  be 

determined. 

The  significant  parameters  and  corresponding  uncertainty  for  the  dilTerent 
configurations  follows  in  Tables  1  through  4.  Graphic  comparisons  of  plate  averaged 
Nusselt  number  versus  Reynolds  number  are  shown  in  Figures  4.1  through  4.7.  These 
graphs  also  include  calculated  uncertainty  ranges  for  the  data,  the  best  fit  correlations 
for  each  configuration,  data  for  comparison  from  other  sources,  and  some  of  the 
correlations  that  have  been  documented  in  research  similiar  to  this  experimentation. 

A.   HEATING  THE  STRAIGHT  UPPER  PLATE 

A  total  of  sixty-nine  data  pairs  were  taken  over  an  approximate  range  of 
Reynolds  numbers  of  6S00  to  20,000  (see  Table  1).  One  means  used  to  analyze  the 
data  was  the  statistics  software,  MINITAB  [Ref  16].  It  was  observed  that  the  data  at 
Reynolds  numbers  below  7000  exhibited  a  relatively  large  deviation  from  the  regression 
in  comparison  to  the  data  for  Reynolds  numbers  above  7000.  This  result  agrees  with 
previous  research  done  by  J.  Wilson  [Ref  17:  p.  50]  which  revealed  that  transition  flow 
in  the  channel  persisted  up  to  a  Reynolds  number  of  approximately  7100;  therefore, 
data  at  Reynolds  numbers  below  7000  was  not  included  in  the  regression  calculated  for 
the  upper  straight  plate. 
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Using  a  least  squares  regression,  a  correlation  for  the  heat  transfer  from  the 
upper  plate  was  determined  to  be: 

N^j  =  0.024Rej^/-^5  {eqn4.1) 

where  Reynolds  and  Nusselt  numbers  are  based  on  the  hydraulic  diameter.  The 
statistical  correlation  of  the  regression  was  0.998.  The  experimental  uncertainty  of  the 
data  was  approximately  ±5%  for  Nusselt  number  and  ±6%  for  Reynolds  number. 

For  comparison,  consider  the  well  known  Dittus-Boelter  correlation  based  on 
hydraulic  diameter  and  a  Prandtl  number  of  0.71: 

Nlij^j  =  0.02Re^jOSO  {eqn4.2) 

At  first  glance  this  small  difierence  between  the  correlations  seems  insignificant; 
however,  a  look,  at  Figure  4.1  shows  that  the  Dittus-Boelter  correlation  overestimates 
the  actual  Nusselt  number  by  about  thirty  to  thirty-seven  percent  over  its  range  of 
applicability  (Re  >  10.000). 

Also  shown  for  comparison  on  Figure  4.1  are  interpolated  values  from  the  Kays 
and  Leung  [Ref  18]  analytical  solution  for  infinite  parallel  plates.  These  values  too  lie 
outside  the  uncertainty  limits  of  the  actual  data,  but  only  overestimate  Nusselt  number 
by  ten  to  fifteen  percent  when  compared  to  the  finite  fiat  plates  of  this  research. 

Data  previously  taken  by  G.  Galyo  [Ref  14:  p.  58]  is  included  in  this  plot.  His 
data  agrees  with  that  of  this  research  to  within  the  bounds  of  uncertainty.  It  must  be 
noted  that  the  data  taken  by  Galyo  at  Reynolds  numbers  above  about  17.000  is  less 
certain  than  the  other  data  due  to  the  difilculty  he  encountered  in  recording  an 
accurate  channel  outlet  temperature  at  high  fiowrates  [Ref  14:  p.  67]. 


B.       HEATING  THE  STRAIGHT  LOWER  PLATE 

A  total  of  sixty-four  data  pairs  were  recorded  over  an  approximate  range  of 
Reynolds  numbers  of  6000  to  19000  (see  Table  2).  As  in  the  previous  case  for  the 
upper  plate,  the  data  recorded  for  Reynolds  numbers  less  than  about  7000  were  not 
included  in  the  least  squares  regression  because  of  their  relatively  high  deviation  from 
the  regression.   The  best  fit  regression  correlation  was: 
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Nu^^  =  O.OlSRej^^O-^s  (eqn4.3) 

The  statistical  correlation  of  the  regression  was  0.998,  and  the  experimental  uncertainty 
of  the  data  was  essentially  the  same  as  that  for  the  upper  flat  plate. 

Figure  4.2  shows  that  over  the  full  range  of  testing  there  is  essentially  no 
difference  in  the  results  for  the  two  flat  plate  cases.  Data  previously  taken  by  G. 
Galyo  [Ref  14:  p.  59]  is  shown  as  well  for  comparison,  and  is  well  within  the 
experimental  uncertainty  of  the  data  recorded  in  this  study  except  at  the  highest 
Reynolds  number  where,  as  mentioned  previously,  G.  Galyo's  data  was  highly 
uncertain. 


C.       HEATING  THE  OUTER  CURVED  (CONCAVE)  PLATE 

Seventy-two  data  pairs  were  taken  over  an  approximate  range  of  Reynolds 
numbers  of  6000  to  20,000  (see  Table  3).  A  least  squares  regression  of  the  data 
resulted  in  the  relation: 

NTi^^^  =  0.05Rej^jO-69  (eqn  4.4) 

The  statistical  correlation  of  the  regression  was  0.999,  and  the  experimental  uncertainty 
was  the  same  as  that  for  the  two  previous  cases. 

Figure  4.3  revealed  that  the  previous  data  reported  by  G.  Galyo  [Ref  14:  p.  60] 
coincides  well  with  the  data  recorded  in  this  study  except,  again,  at  the  highest 
Reynolds  number  where  G.  Galyo  had  diiliculty  determining  duct  outlet  temperature 
accurately.  Also  shown  for  comparison  are  extrapolated  correlations  from  Kreith 
[Ref  12:  p.  1253],  and  Brinich  and  Graham  [Ref  11:  p.  35]. 

In  contrast  to  the  results  for  the  flat  plates  there  was  no  significant  deviation 
from  the  regression  even  at  the  lowest  Reynolds  number  tested  for  the  concave  plate. 
This  supports  previous  research  by  J.  Wilson  [Ref  17:  p.  50]  which  indicated  that  the 
limiting  Reynolds  number  for  transition  flow  in  curved  sections  was  less  than  that  for 
flat  plates  and  approximately  equal  to  5300. 
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The  data  for  the  flat  plates  was  used  as  a  basehne  for  comparison  with  the 
concave  plate,  and  heat  transfer  enhancement  for  the  concave  section  was  reflected  by 
an  increase  in  Nusselt  number  of  about  twenty-one  percent  at  a  Reynolds  number  of 
10,000  to  fourteen  percent  at  a  Reynolds  number  of  20,000  (see  Figure  4.4). 

There  is  a  tendency  to  extrapolate  this  data  to  where  the  correlations  for  the 
concave  plate  and  the  flat  plates  converge  in  order  to  estimate  a  Reynolds  number 
where  the  effect  of  curvature  is  no  longer  significant.  For  these  correlations 
convergence  occurs  at  a  Reynolds  number  of  about  205,000  and  a  Nusselt  number  of 
approximately  232;  however,  research  by  Brinich  and  Graham  indicated  that  the 
enhancement  in  heat  transfer  in  a  concave  section  was  significant  even  at  Reynolds 
numbers  as  high  as  250,000. 

D.       HEATING  THE  INNER  CURVED  (CONVEX)  PLATE 

A  total  of  seventy-three  data  pairs  were  recorded  over  a  range  of  Reynolds 
numbers  of  6000  to  20,000  (see  Table  4).  A  least  squares  fit  of  the  data  gave  the 
relation: 

^u^^  =  O.OlTReJ-'^^  (eqn4.5) 

at  a  statistical  correlation  of  0.997,  and  an  experimental  uncertainty  the  same  as  for  the 
previous  cases. 

Figure  4.5  revealed  that,  as  was  the  case  for  the  concave  plate,  there  was  no 
significant  deviation  from  the  regression  even  at  the  lowest  Reynolds  number  tested.  A 
comparison  of  this  data  with  that  reported  by  G.  Galyo  [Ref  14:  p.  61]  showed  that 
there  is  agreement  to  within  the  bounds  of  experimental  uncertainty  except  at  the 
highest  Reynolds  number.  For  comparison,  data  reported  by  Kreith  [Ref  12:  p.  1253], 
and  Brinich  and  Graham  [Ref  11:  p.  36]  for  a  convex  surface  were  included  in  the 
graph. 

Using  the  flat  plate  data  as  a  reference  of  comparison,  the  reduction  in  heat 
transfer  on  the  convex  plate  was  manifest  in  a  Nusselt  number  consistently  six  percent 
lower  than  that  of  the  flat  plates  over  the  entire  range  of  How  rates.  This  did  not 
represent  a  large  difference  between  the  convex  and  flat  plate  cases;  however,   six 
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percent  was  at  least  one  percent  more  than  the  conservatively  calculated  experimental 
uncertainty  of  the  data.  Therefore,  this  data  implied  that  a  decrease  in  heat  transfer  on 
a  convex  surface  should  be  expected,  though  probably  small. 

A  comparison  of  the  convex  and  concave  surfaces  (see  Figure  4.6)  pointed  out 
that  the  concave  surface  heat  transfer  is  approximately  thirty  percent  greater  than  that 
for  the  convex  surface  at  a  Reynolds  number  of  8000,  and  about  twenty  percent  greater 
at  a  Reynolds  number  of  20,000. 

E.       EFFECT  OF  BUOYANCY 

Tables  3  and  4  include  data  taken  when  the  test  channel  was  oriented  such  that 
free  convection  aided  forced  convection.  Graphical  comparison  for  this  case  is  shown 
in  Figure  4.7. 

Data  taken  for  this  configuration  was  within  the  bounds  of  experimental 
uncertainty,  and  showed  no  pattern  in  difiering  from  the  case  where  the  direction  of 
free  convection  opposed  that  of  forced  convection.  These  results  indicated  that  there 
was  no  significant  buoyant  effect  over  the  full  range  of  tlow  rates  studied. 
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TABLE  1 
SUMMARY  OF  STRAIGHT  UPPER  WALL  RESULTS 


R^hd 

^\6 

^^^hd^^hd 

6841 

17.30 

±0.0560 

7856 

19.92 

8588 

21.46 

9327 

23.06 

10168 

24.75 

±0.0576 

10936 

26.18 

11947 

28.11 

12082 

28.38 

12981 

29.25 

13969 

30.96 

±0.0571 

14997 

32.66 

15950 

34.62 

17169 

36.85 

17839 

37.30 

±0.0570 

18927 

39.01 

20008 

40.27 

±0.0570 

^^\<i^\6 


±0.0481 


±0.0467 


±0.0485 


±0.0495 
±0.0509 
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TABLE  2 
SUMMARY  OF  STRAIGHT  LOWER  WALL  RESULTS 


Rsd 

^■^H.d 

SR^hd^^hd 

^^■^d^'^d 

6071 

15.03 

±0.0581 

±0.0452 

7065 

17.85 

8116 

19.72 

±0.0573 

±0.0458 

8968 

21.26 

9997 

23.63 

10934 

25.35 

±0.0571 

±0.0464 

12014 

27.17 

12984 

29.11 

14047 

30.89 

±0.0571 

±0.0473 

14989 

32.67 

15981 

34.10 

16984 

35.45 

±0.0570 

±0.0479 

18035 

36.21 

19002 

37.46 

19458 

39.83 

±0.0570 

±0.0474 
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TABLE  3 
SUMMARY  OF  CONCAVE  WALL  RESULTS 


R^hd 

^^^hd 

SR^hd'I^^hd 

^^^^d'^'^d 

5884 

20.04 

±0.0582 

±0.0448 

6025* 

19.98* 

7050 

22.89 

8016 

24.73 

±0.0573 

±0.0450 

8044* 

24.80* 

9094 

26.89 

10025* 

28.08* 

10045 

28.71 

11023 

30.72 

±0.0571 

±0.0451 

11977 

32.72 

13028 

34.90 

13993 

36.19 

±0.0571 

±0.0459 

14982* 

38.05* 

15021 

37.95 

15996 

40.11 

17053 

42.51 

±0.0570 

±0.0464 

17943 

42.90 

19169 

45.26 

20189 

46.80 

±0.0571 

±0.0478 

*   Indicates  measurements  that  were  taken  with  the  test  channel  inverted 
such  that  the  directions  of  free  and  forced  convection  were  coincident. 
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TABLE  4 
SUMMARY  OF  CONVEX  WALL  RESULTS 


R^hd 

^Ad 

SR^hd^^hd 

^^Ad^^'^d 

6065 

15.04 

±0.0581 

±0.0492 

6088* 

15.15* 

7050 

16.85 

8009 

18.87 

±0.0574 

±0.0498 

8034* 

18.88* 

9058 

19.90 

10020 

22.12 

10040* 

22.48* 

11054 

24.32 

±0.0571 

±0.0525 

12003 

24.92 

12988 

27.83 

13944 

29.42 

±0.0571 

±0.0545 

14871* 

30.57* 

15003 

30.72 

15981 

32.20 

17009 

32.84 

±0.0570 

±0.0517 

17898 

35.43 

19015 

36.25 

20210 

37.73 

±0.0571 

±0.0524 

*   Indicates  measurements  that  were  taken  with  the  test  cliannel  inverted 
such  that  the  directions  of  free  and  forced  convection  were  coincident. 
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Figure  4.1     Straight  Upper  Plate  Results. 
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Figure  4.2     Straight  Lower  Plate  Results. 
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Figure  4.3     Concave  Plate  Results. 
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Figure  4.4    Comparison  of  Straight  Upper  and  Concave  Wall  Results. 
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Figure  4.5     Convex  Plate  Results. 
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Figure  4.6    Comparison  of  Concave,  Convex  and  Straight  Plate  Data. 
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Figure  4.7    Convex  Concave  Buoyancy  Comparison. 
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V.  DISCUSSION 

A.       FLAT  PLATES 

The  experimental  results  for  the  upper  and  lower  straight  sections  of  the  channel 
were  essentially  the  same,  and  are  not  considered  separately  in  this  discussion. 

The  Dittus-Boelter  correlation  for  Xusselt  number  based  on  the  hydraulic 
diameter  of  the  test  channel  and  a  Prandtl  number  of  0.71  yielded  Nusselt  numbers 
thirty  to  thirty-seven  percent  higher  than  those  recorded  in  this  experimentation.  The 
Dittus-Boelter  correlation  is  an  empirical  approximation  that  was  developed  to 
estimate  local  Nusselt  numbers  for  circular  tubes.  The  correlation  itself  is  only  an 
approximation,  and  when  modified  for  non-circular  ducts  where  convection  coefficients 
van.'  about  the  peripher\',  the  approximation  becomes  less  certain.  Here  the  correlation 
was  used  to  estimate  the  average  Nusselt  number  over  the  length  of  the  heated  section. 
To  a  first  approximation  the  Dittus-Boelter  correlation  can  be  used  to  estimate  the 
average  Nusselt  number  because  the  ratio  of  the  plate  length  to  the  hydraulic  radius 
was  greater  than  ten:  however,  this  assumption  adds  to  the  uncertainty  of  the  resulting 
convection  coefficient. 

The  solution  for  parrallel  plates  of  infinite  extent  developed  by  Kays  and  Leung 
gave  Nusselt  numbers  closer  to  the  experimental  results,  but  was  still  high  by  ten  to 
fifteen  percent,  though  ten  to  fifteen  percent  accuracy  when  estimating  heat  transfer 
coefficients  in  the  turbulent  flow  regime  is  not  unusual.  Even  though  the  straight 
section  of  the  test  channel  approximated  flat  parallel  plates  of  infinite  extent,  some 
sidewall  influence  on  flow  was  expected  and  may  account  for  a  part  of  the  disparity  of 
the  comparison.  The  graphical  con^parison  shown  in  Figure  4.1  was  developed  by 
linear  interpolation  of  the  data  from  [Ref  18:  p.  552],  though  there  was  no  evidence  to 
support  the  assumption  that  the  results  of  the  Kays  and  Leung  solution  are  a  linear 
function  of  either  Prandtl  or  Reynolds  number;  therefore,  additional  inaccuracy  may 
have  been  introduced  into  the  comparison. 
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B.       CURVED  PLATES 

Experimental  and  analytical  investigation  of  flow  and  heat  transfer  on  concave 
and  convex  surfaces  in  a  rectangular  duct  of  high  aspect  ratio  has  thus  far  received 
little  attention;  consequently,  there  is  scant  previous  work  available  for  comparison. 
Some  results  are  available  for  channels  of  modest  aspect  ratios,  but  these  do  not 
generally  compare  well  due  to  the  significant  sidewall  effect  introduced  at  low  aspect 
ratios. 

Studies  of  the  flow  and  heat  transfer  characteristics  in  high  aspect  ratio  channels 
has  been  done  at  the  Naval  Postgraduate  School  [Refs.  9,14,17,19,20,21,22].  Results  of 
this  previous  research  were  comparable  to  the  results  of  this  experimentation  and  were 
highly  depended  upon  to  validate  these  results.  Only  Galyo's  work,  is  shown  for 
comparison  in  this  report,  but  his  results  agreed  within  the  limits  of  experimental 
uncertainty  with  those  of  Wilson  [Ref  17],  and  Daughety  [Ref  22].  The  major 
improvements  of  this  experimentation  over  that  previously  done  were  the  decreased 
experimental  uncertainty,  especially  at  the  highest  Reynolds  number,  and  the 
preponderance  of  data  that  permitted  empirical  correlations  to  be  determined  to  a  very 
high  level  of  confidence. 

In  1955  Kreith  [Ref  12:  p.  1253]  investigated  the  influence  of  curvature  on  heat 
transfer  in  a  rectangular  channel  of  aspect  ratio  equal  to  three.  For  a  Prandtl  number 
of  0.71,  graphical  comparison  with  the  results  recorded  in  this  experimentation  showed 
that  Kreith's  results  were  approximately  thirty  percent  lower  (see  Figures  4.3  and  4.5); 
however,  Kreith's  test  apparatus  not  only  had  a  much  lower  aspect  ratio,  it  was  also 
heated  by  a  heating  element  that  was  centered  between  the  upper  and  lower  walls  of 
the  channel  instead  of  heating  the  walls  themselves. 

Research  done  in  1977  with  a  channel  of  aspect  ratio  six  by  Brinich  and  Graham 
[Ref  11:  pp.  35-36]  was  shown  for  comparison  in  Figures  4.3  and  4.5.  The  range  of 
Reynolds  numbers  over  which  they  performed  their  experiments  was  approximately 
84,000  to  253,000.  For  the  concave  surface,  Brinich's  and  Graham's  data  extrapolated 
to  a  Reynolds  number  of  15,000  yielded  a  Nusselt  number  about  fifty  percent  higher 
than  that  recorded  in  this  experiment.  Likewise,  an  extrapolation  to  a  Reynolds 
number  of  15.000  for  the  convex  plate  showed  that  there  was  about  a  thirty  percent 
difference  in  the  results  of  the  two  experiments.    These  difTerences  are  rather  large,  but, 
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again,  the  two  studies  diflered  greatly  in  that  the  aspect  ratios  and  the  range  of 
Reynolds  numbers  investigated  were  much  difierent. 

Similiar  research  was  done  in  1978  by  L'chida  and  Koizumi  [Ref.  13:  pp. 
1714-1716],  with  channels  of  aspect  ratio  up  to  four;  however,  the  full  report  has  not 
been  translated  into  English.  From  figures,  graphs,  and  equations  included  in  their 
report  it  appears  that  their  test  channel  was  heated  along  its  entire  periphery'  and  that 
the  published  results  were  in  the  form  of  an  average  Nusselt  number  for  the 
combination  of  a  convex  and  concave  wall. 

It  is  unlikely  that  an  average  heat  transfer  coefficient  or  Nusselt  number  can  be 
determined  when  both  the  convex  and  concave  walls  are  heated  by  simple 
superposition  of  the  individual  correlations  for  the  convex  and  concave  walls  because 
of  the  complex  nature  of  the  fiow  in  a  curved  channel.  Therefore,  there  was  no 
manner  in  which  the  results  detailed  in  Ref  13  could  justifiably  be  compared  to  the 
results  of  this  experimentation. 

Even  though  the  test  equipment  was  such  that  a  pair  of  heaters  could  be 
energized  to  supply  a  constant  heat  flux  on  opposite  walls  of  the  channel,  no 
conclusive  experiments  could  be  performed  to  determine  an  average  Nusselt  number 
for  this  configuration. 

Referring  to  (eqn.  3.3)  and  denoting  hj  as  the  average  heat  transfer  coefficient  for 
single  plate  heating,  and  h2  as  the  average  heat  transfer  coefficient  when  heating  a  pair 
of  opposite  walls,  then: 

h,         [(T      +  T  .)  2  -  T. ,.  ] 

^  ^  iLwo uv^ bikJ  (eqn  5.1) 

^2  (T^,,  -  T,„) 

where  T  ,  represents  the  average  temperature  when  one  wall  (inner  or  outer)  is  heated, 
and  the  heat  flux  from  both  the  inner  and  outer  plates  was  assumed  equal  when 
heating  a  pair  of  plates.  The  average  wall  temperature  term  in  the  numerator  of  (eqn. 
5.1)  can  be  recast  as: 

(T      +  T  .)  2  =  T  .    4-  5T  (eqn  5.2) 

^    wo  vvr  vvn  ^    ^ 
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where  6T  represents  the  diflerence  between  the  average  temperature  when  heating  a 
pair  of  plates  and  the  single  plate  average  temperature.  Combining  (eqn.  5.1)  and  (eqn. 
5.2)  yields: 

J_  =^^vh bik (        53) 

_  Ij^illLWbikL 


1 

I  +  6T  AT 

i 

(Eqn.  5.3)  reveals  that  only  when  the  temperatures  of  both  plates  are  the  same,  such 
that  6T  =  0,  when  heating  a  pair  of  plates,  can  a  true  average  heat  transfer  coeHicient 
for  both  walls  be  obtained.  Strictly  speaking  the  temperature  at  all  corresponding 
points  opposite  each  other  on  the  plates  should  be  equal  to  yield  an  average  heat 
transfer  coelTicient  for  the  combination  of  both  walls.  The  heaters  used  in  this 
experimentation  were  constant  heat  flux  heaters,  and  temperatures  at  corresponding 
points  on  opposite  walls  could  not  be  matched;  therefore,  conclusive  data  for  the  case 
of  heating  a  pair  of  plates  could  not  be  gathered  using  the  existing  equipment. 

Though  local  heat  transfer  coefficients  and  Nusselt  numbers  were  calculated  at 
the  leading  edge,  at  mid-plate,  and  at  the  trailing  edge  of  the  heated  plates,  no 
conclusive  pattern  was  evident  in  the  data.  The  same  was  true  for  local  temperatures. 
More  conclusive  results  could  not  be  obtained  because  there  were  too  few 
thermocouples  spaced  too  far  apart;  that  is.  the  spanwise  thermocouple  spacing  was 
greater  than  the  diameter  of  the  Taylor-Gortler  vortices  that  formed. 
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VI.  CONCLUSION 

The  results  of  this  experimentation,  detailed  in  Chapter  IV,  indicated  that 
secondare'  flow  in  a  curved  rectangular  channel  in  the  form  of  Taylor-Gortler  vortices 
influenced  the  heat  transfer  on  the  concave  and  convex  surfaces  over  a  range  of  flow 
corresponding  to  Reynolds  numbers  from  about  6000  to  20000. 

Correlations  of  Nusselt  number  as  a  function  of  Reynolds  number  for  the  two 
flat  plates,  the  convex  plate,  and  the  concave  plate  were  calculated  to  a  high  level  of 
certainty  and  showed  that  heat  transfer  was  enhanced  by  approximately  twenty  percent 
on  the  concave  surface  as  compared  to  the  flat  surfaces.  Similiarly,  heat  transfer  over 
the  convex  surface  was  reduced  by  approximately  six  percent  in  comparison  to  the  flat 
surfaces.  Furthermore,  the  enhancement  or  reduction  in  heat  transfer  on  the  curved 
plates  was  noted  to  be  greatest  at  the  lowest  Reynolds  number  tested,  and  tapered  ofl" 
with  increased  flow  rate. 

Comparisons  made  of  the  results  obtained  in  this  experimentation  to  that  of 
other  pertinent  research  supported  the  conclusion  that  heat  transfer  is  enhanced  on  a 
concave  surface.  Experiments  performed  for  the  purpose  of  determining  the  relative 
significance  of  buoyancy  induced  flow  indicated  that  even  at  the  lowest  Reynolds 
number  of  6000.  buoyant  free  convection  was  negligible  compared  to  forced 
convection. 


46 


VII.  RECOMMENDATIONS 

A.       EQUIPMENT 

1.  Small  oscillations  in  ambient  temperature  (±1°)  noticeably  afTected 
temperatures  in  the  channel  and  contributed  to  the  difTerences  in  the  results  in 
successive  data  taking  runs.  Maintaining  close  control  over  ambient 
temperature  would  make  experimental  results  more  reproducible  and  decrease 
the  time  necessary  to  gather  data. 

2.  Existing  equipment  allowed  for  testing  only  up  to  a  Reynolds  number  of  about 
20000,  though  the  corresponding  Reynolds  numbers  for  flow  through  a  gas 
turbine  are  much  higher.  A  larger  test  channel  and  higher  capacity  air 
compressor  would  be  required  for  experimentation  at  higher  flow  rates. 

3.  A  great  amount  of  time  was  spent  waiting  approximately  four  hours  for  the  test 
channel  to  reach  steady  state  after  first  energizing  the  heater.  Installing  guard 
heaters  would  significantly  reduce  the  warm  up  time  to  steady  state. 
Additionally,  guard  heaters  would  prevent  heat  loss  through  the  insulation 
surrounding  the  heaters. 

4.  Pressure  sensors  that  could  be  read  directly  by  the  automatic  data  acquisition 
system  several  times  during  the  course  of  a  data  run  would  increase  the 
accuracy  of  flow  rate  measurements  and  permit  successive  data  runs  to  be  taken 
without  the  presence  of  an  operator. 

5.  More  thermocouples  should  be  installed  closer  together  to  improve  the  accuracy 
of  temperature  measurements  and  to  facilitate  studying  the  pattern  of  local 
variations  in  temperature  and  heat  transfer  in  the  spanwise  and  flow-wise 
directions. 

6.  The  single  most  significant  contributor  to  experimental  uncertainty  in  the 
calculated  average  Nusselt  number  and  hydraulic  Reynolds  number  was  the 
uncertainty  in  the  channel  cross  section.  The  uncertainty  in  experimental 
results  could  be  greatly  reduced  if  a  channel  were  constructed  paying  special 
attention  to  the  cross  section  dimensions. 
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B.       FURTHER  RESEARCH 

1.  Experimentation  similiar  to  this  but  at  higher  flow  rates  should  be  done  so  that 
a  large  data  base  could  be  established  and  heat  transfer  correlations  could  be 
calculated  for  a  wide  range  of  application. 

2.  Experiments  at  lower  flow  rates  where  buoyancy  is  significant  should  be  done 
so  that  a  suitable  length  parameter  can  be  deduced  for  calculating  the  Grashof 
number  for  curved  channel  flow. 

3.  Different  curvatures  should  be  studied  to  determine  the  effect  of  the  degree  of 
curvature  on  heat  transfer  in  a  curved  channel. 

4.  Thermostatic  control  should  be  included  in  the  experimental  equipment  that 
could  maintain  equal  plate  average  temperatures  when  heating  a  pair  of  heaters 
on  opposite  walls  so  that  experimental  results  for  channel  average  heat  transfer 
coefficients  could  be  gathered. 
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APPENDIX  A 
DATA  REDUCTION  PROGRAM 


1000  !   FILE  NAME:  MAIN 

1010  !   LAST  REVISION  DATE:   19  OCT  1986 

1020  !   EDITED  BY:   LT.  J.R.  HAWK   III 

1030! 

lOAO  !   PROGRAM  FOR  GATHERING  AND  REDUCING  DATA  WHEN  HEATING  ONE  COPPER  PLATE, 

10S0  !   WITH  S  MINUTES  BETWEEN  RUNS.   PROGRAM  IS  MODIFIED  TO  ACCOUNT  FOR 

lOGO  !   BROKEN  THERMOCOUPLE  CHANNELS  5,  16,  21,  AND  26  WHICH  EFFECT  OUTER 

1070  !   CURVED  HEATER.  THERMOCOUPLE  CHANNEL  47  WtllCH  EFFECTS  INNER  CURVED 

1080  !   HEATER,  AND  THERMOCOUPLE  CHANNEL  11  WHICH  EFFECTS  OUTER  STRAIGHT  HEATER 

1090! 

1100  COM  /Co/  D(7) ,Aa(76,2) 

1110  DIM  Emf (79) ,T(79) 

1120! 

1130  I   CORRELATION  FACTORS  FOR  EMF  TO  DEGREES  CELSIUS  FOR  CU-CO  T/C. 

1  W)0! 

11&0  DATA  0.  10086091 ,25727.94369, -76734S.829S.7802559S. 81 

1160  DATA    -9247486S89,6.97688E  +  11,-2.66192E  +  13,3.94078EH4 

1170  READ   D(») 

1180  PRINIER  IS  701 

1190  BFFP 

1200  INPUT  "ENTER  MONTH,  DAY,  AND  TIME  ( MM:DD:HH:MM:SS) " , T ime$ 

1210  OUTPUT  709;"TD";Time$ 

1220  BEEP 

1230  INPUT  "ENTER  INPUT  MODE  ( 0=3054A. 1 -FILE )", Ira 

1240  IF  Im-1  THEN 

1250  BEEP 

1260  INPUT  "ENTER  NAME  OF  EXISTING  DATA  FILE" ,01df 1 le$ 

1270  PRINT  USING  "20X, ""These  Results  Are  From  Data  File:   "  "  ,  1 OA" ;01df i le$ 

1280  PRINT 

1290  PRINT 

1300  ELSE 

1310  BEEP 

1320  INPUT  "ENTER  NAME  OF  NEW  DATA  FILE?" ,Neuf i le$ 

1330  PRINT  USING  "20X,""Data  For  This  Run  Is  Stored  In  Data  File:   " "  .  1 0A"  ;Neujf 

ileS 

1340  PRINT 

1350  PRINT 

1360  END  IF 

1370! 

1380  !   ENTER  THE  CHANNEL  CONFIGURATION. 

1390! 

HOG  PRINT  USING  "9X,""THE  FOLLOWING  CHANNEL  CONFIGURATION  WAS  SELECTED:""" 

1410  PRINT 

1420  BEEP 

1430  INPUT  "SELECT  HEATER  TYPE  (0  =  STRAIGHT,  1  =CURVED ) " . I  type 

1440  IF  ltype  =  0  THEN  PRINT  USING  "  1  2X . " "Hea 1 1 ng  STRAIGHT  Test  Section.""" 

1450  IF  Itype=1  THEN  PRINT  USING  "  1 2X , " "Heat i ng  CURVED  Test  Section.""" 

1460  BEEP 

M70  INPUT  "SELECT  HEATER  POSITION  (0-OUTER,  1  =  INNER )", Ipos 

1480  IF  Ipos-0  THEN  PRINT  USING  "  1  2X , " "Hea 1 1 ng  OUTER  Plate. 

1490  IF  Ipos-I  THEN  PRINT  USING  "  12X , " "Heat i ng  INNER  Plate.""" 

1500  PRINT 

1S10! 
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,.0   .   ENTERS  riLES  UUH  UC   C«UB.«T:ON  COEEEICIENTS  B«EO  ON  H«TE.  T.PE. 

lis;  r,  1-:^  ]T.  ^^  lEiifc  IS  :irv-isir. 

j^lS      SEslG^SFUec   TO   . 

!6?3'    I      INSERT   DflTfl  STORAGE  DISK. 

\lf,      ['Arizs'"  ^e.Ule^,^0 

1690  END  IF 

1700   ENTER  709:Time$ 

1?^„°   ^^^S?aO-20X.--Non.^.  Oav.  an.  T.»e:",.--.T.»eS 

1730   PRINT 

l^sS'-  !   FOR  DAT.  RUN  READ  AND  STORE  ALL  RAU  E.F  VALUES. 

y}l  JUSt'/Ss'-AR  AFOO  AL79- 

90  OUTPUT  722;;F1  R1  T1  Zl  FLl 

leOO  FOR  I'O  TO  9 

1810  OUTPUT  709;  AS 

1820  ENTER  722;Emt(I) 

\lll      SuTPuf  L*FiIe;En,f(*) 

18S0  ELSE 

I860  ENTER  (^F  i  le'.Emf  (  *  » 

l^ir  PRINT  USING  -SX.-THE  FOaONING  DATA  UERE  RECORDED: 

1830   PRINT 

iiir;lHEri°?3rj^^i;:^.s^w.THEN 

1890   PRINT  USING   l^X.   Heater   o  ^  ^^    ^^  ^^  __  ^^^^ ;Emt(39>Ia 

.^SoO   PRINT  USING  ■•l2X.-"Res.stor  Voltage  .Vr)     '    •-  •   '  ^   ^^^^^^^  ^^^^ 

^Sin   IF  IPOS.)  THEN  PRINT  USING  ••  1 2K , '• -P  rec  .s.on   es.s  to.  ^   ^  ^^^^^^  ^^^^ 
f^:;;-  if  Ipos^  then  print  using  ■.2X.-'Prec.s.onRes.stor  .Rpc 

s ) 

2030  PRINT 
20AO  END  IF 

11'  !  'if  T/CQUT  of  range  PRINT  A  QUE. 
2080! 


50 


2090 

Aemt - 

Etnf  (I) 

2100 

IF  Aemf<7.09E-4  OR  Aemf >3. 630E- 

-3  THEN 

2110 

T(I)  = 

999.99 

2120 

ELSE 

2130 

CALL 

Tv/sv(  Emi 

r(I).Tt) 

2M0 

IF  I- 

79  THEN 

2150 

T(I)- 

Tt 

2160 

GOTO 

2210 

2170 

END  IF 

2180 

T(I)- 

FNTcorr 

(Tt,I) 

2190 

END  IF 

2200 

END  IF 

2210 

NEXT 

I 

2220! 

2230 

!   PRINT  ALL 

T/C  READINGS. 

2240! 

2250 

PRINT 

USING  ' 

•12X. ""OUTER  PLATE; 

,  M  !•  M 

2260 
6 "  - 
2270 

PRINT 

USING  ' 

"14X.""T/C  Number: 

1 

2 

PRINT 

USING  ■ 

•14X,""Temp  (C): 

"".6 (3D, 

,DD.2X) 

(5) 

2280 

«  Til  n  II 

PRINT 

USING  ' 

•MX.'-T/C  Number: 

7 

8 

1  C 

2290 

PRINT 

USING  ' 

'MX,  ""Temp  (C): 

"".6(3D. 

,DD.2X) 

T(1  n 

2300 

PRINT 

USING  ' 

'14X,""T/C  Number: 

13 

14 

18 

2310 

PRINT 

USING  • 

'14X,""Temp  (C): 

"",6 (3D. 

,DD,2X) 

16) .T(17) 

2320 

PRINT 

USING  ' 

'14X,""T/C  Number: 

19 

20 

24 

2330 

PRINT 

USING  ' 

"MX,  ""Temp  (C): 

"-.6 (3D. 

DD,2X) 

22),T(23) 

2340 

PRINT 

USING  ' 

'14X,""T/C  Number: 

25 

26 

3C""" 

2350 

PRINT 

USING  " 

'14X,""Temp  (C): 

"-.5(30. 

DD,2X) 

28).T(29) 

2360 

PRINT 

2370! 

2380 

!   CALCULATE 

AVERAGE  TEMPERATURE  OF  OUTER  HEAT 

2390! 

2400 

Sum  =  0 

^ 

2410 

IF  It' 

/pe-0  THEN 

2420 

FOR  \ 

=0  TO  2S 

) 

2421 

IF  1  = 

11  THEN 

2440 

2430 

Sum  =  Si 

jm+T(I) 

2440 

NEXT 

I 

2450 

Tu)o  =  Sum/29 

24bO! 

2470 

!   CORRECTION 

1  FOR  BROKEN  T/C  CHANNELS  5. 

16,  2 

ATER. 

2480! 

2490 

ELSE 

2500 

FOR  1  = 

=0  TO  2S 

) 

2510 

IF  I  =  i 

S  OR  1=16  OR  1-21  OR  1=26 

THEN  2530 

2520 

Sum  =  Si 

jm+T( I) 

2530 

NEXT  ; 

I 

2540 

Tijo  =  Si 

jm/25 

2550 

END  IF 

2560 

PRINT 

USING  " 

12X, ""OUTER  INSULATION:""" 

3       4       5 
:T(0) ,T( 1 ),T(2) ,T(3),T(4),T 

9      10      11 
:T(6),T(7),T(8),T(9) ,T(10), 

15      16      17 
T(12) ,T(13) ,T( 14) ,T(15) ,T( 

21      22      23 
T( 18) .T( 19),T(20),T(21 ) ,T{ 

27      28      29 
T(24).T(25),T(26),T(27),T( 


AND  26  IN  OUTER  CURVED  HE 
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2570   PRINT  USING  -MX,""T/C  Number:     31  (Tinsol)    32  (Tinso2)    33  (Tinso3)" 

2B80  PRINT  USING  "MX,""Tenp  (C):        " " ,  3(  3D.DD,8X ) " ;  T(  30  )  ,  T(  31  )  ,  T(  32) 

2S90   Tinsol-T(30) 

2600   Tinso2-T(3n 

2610   Tinso3-T(32) 

2620   PRINT 

2630  PRINT  USING  " 12X, ""ORIFICE  TEMP  (Torf)  •  "",3D.DD.-"  (C) '"■ " ; 1(33) 

26-^0   PRINT 

2550  PRINT  USING  "12X, ""INLET  TEMPERATURE:""" 

2660   PRINT  USING  "1^X,""T/C  Number:     35      36      37      38""" 

2670   PRINT  USING  "MX. ""Temp  (C):     " "  ,  4(  3D  .  DD  .  2X  )  ' ;  T(  34  )  ,  T(  35  )  .  T(  36  )  ,  T(  37  ) 

2680   PRINT 

2690! 

2700   !   CALCULATE  AVERAGE  INLET  TEMPERATURE. 

2710! 

2720   Sum-0. 

2730   FOR  I-3A  TO  37 

27^0   Sum=Sum+T(I) 

2750   NEXT  I 

2760   Tin-Sum/4 

2770  PRINT  USING  "12X, ""INNER  PLATE:""' 

2780   PRINT  USING  "14X.""T/C  Number:     41      42      43      44      45 

46 

2790   PRINT  USING  "MX. ""Temp  (C):     "  "  ,6(  3D.  DD,2X  )  " ;  T(  40  )  .  T(  4  1)  ,  T(  42  )  ,  T(  43)  ,  T( 

44)  .T(45) 

2800   PRINT  USING  "14X,""T/C  Number:     47      48      49      50      51 

52 

2810  PRINT  USING  "14X,""Temp  (C):     " " ,6 ( 3D . DD,2X ) " ; T( 46) . T( 47 ) . T( 48 ) , T( 49) , T( 

50) .T(51 ) 

2820   PRINT  USING  "14X,""T/C  Number:     53      54      55      56      57 

53 

2830  PRINT  USING  "I4X.""Temp  (C):     " " ,S( 3D. DD.2X ) " ; T (52 ) , T( 53 ) . T( 54 ) . T(55) . T ( 

5b).T(57) 

2840  PRINT  USING  "14X,""T/C  Number:    59      60      61      62      53 

64 

2850  PRINT  USING  "14X.""Temp  (C):     " " . 6 ( 3D. DD, 2X ) " : T(58 ) . T ( 59 ) . T( 60 ) . T ( 61 ) , T( 

6"' )  T  (  6 3  ) 

2860   PRINT  USING  "14X.""T/C  Number:     65      66      57      68      69 

70 

2870   PRINT  USING  "14X.""Temp  (C):     " " .6( 3D. DD,2X ) " ; T ( 64 ) . T( 65 ) . T(66 ) , T ( 67 ) , T ( 

68)  .T(6'3) 

2880   PRINT 

2890! 

2900   !   CALCULATE  AVERAGE  TEMPERATURE  OF  INNER  HEATED  WALL. 

2910! 

2920   Sum-0. 

2930   IF  ltype=0  THEM 

2940   FOR  1-40  TO  69 

2950   Sum=Sum+T(I) 

2960   NEXT  I 

2970   Tiji=Sum/30 

2980! 

2990   !   CORRECTION  FOR  BROKEN  T/C  CHANNEL  47  IN  INNER  CURVED  HEATER. 

3000! 

3010  ELSE 

3020  FOR  1=40  TO  59 

3030  IF  1=47  THEN  3050 

3040  Sum=Sum+T( I) 

3050  NEXT  I 
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30GO 
3070 
3080 
3090 

II  H 

3100 

3110 

3120 

3130 

31^0 

3150 

3160 

3170 

3180 

3190! 

3200 

3210! 

3220 

3230 

3240 

3250 

3260 

3270! 

3280 

3290! 

3300 

3310! 

3320 

3330! 

3340 

3350! 

3360 

3370! 

3380 

3390 

3400 

DD."" 

3410 

DD,"" 

3420 

DD."" 

3430 

DD . " " 

3440' 

DD."" 

3450 

DD."" 

3460! 

3470 

34:j0! 

3490 

3500 

3510 

3520 

3530 

DD."" 

3540! 

3550 

3560! 

3570 


Tui=Sum/29 

END  IF 
PRINT  USING  "12X. ""INNER  INSULATION:""" 
PRINT  USING  "14X,""T/C  Number:     71  (Tinsil) 


72  (Tinsi2)    73  (Tinsi3) 

PRINT  USING  "MX, ""Temp  (C):        ""  .3(3D.DD,8X) " ;  T(  70  )  ,  T(  71  )  .  T(  72) 

Tinsi1-T(70) 

Tinsi2"T(71  ) 

Tinsi3-T(72) 

PRINT 

PRINT  USING  "12X. ""OUTLET  TEMPERATURE:""" 

PRINT  USING  "14X,""T/C  Number:     74      75      7S      77""" 

PRINT  USING  "14X,""Temp  (C):     " " , 4( 3D. DD,2X) " ; T( 73 ) , T( 74 ) . T ( 75 ) , T( 76 ) 

PRINT  USING  "*.*" 

!   CALCULATE  AVERAGE  OUTLET  TEMPERATURE. 

Sum-0. 

FOR  1-73  TO  76 

Sum"Sum+T( I ) 

NEXT  I 

Tou  t-Sum/4 

!   CALCULATE  DIFFERENCE  BETUEEN  CHANNEL  INLET  AND  OUTLET  TEMPERATURES. 

Tdif f-Tout-Tin 

!   CALCULATE  CHANNEL  BULK  TEMPERATURE. 

Tblk-(Tin+Tout)/2 

!   PRINT  CALCULATED  TEMPERATURES. 


PRINT  USING  "9X,""THE  FOLLOUING  TEMPERATURES  WERE  CALCULATED: 

PRINT 

PRINT  USING  "12X. ""Average  Outer  Wail  Temperature  (Two) 


(C 

PRINT 
( C) " " 
PRINT 
(  C  )  " " 
PRINT 
(  C  )  " " 
PRINT- 
(  C)  " " 
PRINT 
(  C  ) " " 


" ;  Tuio 

__USING    "12X. 
" ;  T  ij  1 

USING    "12X, 
" ;  Tout 

USING    "12X. 
"  ;  T  1  n 

USING    "12X, 
";Tdiff 

USING    "12X. 
";Tblk 


""Average  Inner  Wall  Temperature  (Tuii) 

""Average  Outlet  Temperature  (Tout) 

""Average  Inlet  Temperature  (Tin) 

""Channel  Inlet  and  Outlet  Temp  Diff  (Tdiff) 

""Average  Bulk  Temperature  (Tblk) 


,3D, 
,3D. 
,3D, 
,3D. 
.3D. 
,3D. 


!   CALCULATE  DELTA  T  (Tdel)  FOR  USE  IN  CALCULATING  H  LATER, 

7delo=Tuo-Tbik 

Tdeii  =  Tiji-Tblk 

IF  lpos=0  THEN  Tdel'Tdelo 

IF  Ipos=l  THEN  Tdel=Tdeli 

PRINT  USING  "12X.""Mean  Temperature  Difference  (Tdel) 

(C) ;Tdei 

!   CALCULATE  DIFFERENCE  BETWEEN  OUTER  AND  INNER  WALL. 

Tudiff-ABS( Two-Tui) 


,3D, 
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3S80  PRINT  USING  " 12X,  ""Outer  and  Inner  Wall  Temp  Difference  (Tudiff)   -  "",3D. 

DD,""  (C) ;Tuidiff 

3590! 

3600  !   CALCULATE  LOCAL  CHANNEL  TEMPERATURES  AT  INLET.  MIDDLE  AND  OUTLET. 

3610! 

3520  Tcin-Tin+( .5*Tdiff/12) 

3630  Tcmid-Tin+(4.5*Tdiff/12) 

35-^0  lcout  =  Tin+(  11  .5*Tdiff/12) 

3650  PRINT  USING  " 12X ." "Average  Local  Channel  Temp  at  Inlet  (Tcin)     -  "".30. 

DD.""  (C)""";Tcin 

3660  PRINT  USING  " 12X , ""Average  Local  Channel  Temp  at  Middle  (Tcmid)    -  ",3D. 

DO.""  (C)""":Tcmid 

3670  PRINT  USING  " 12X, ""Averaae  Local  Channel  Temp  at  Outlet  (Tcout)    -  "",30. 

DD,""  (C) ; Tcout 

3680! 

3690  !   CALCULATE  AVERAGE  LOCAL  HEATED  WALL  TEMPERATURE  AT  INLET. 

3700! 

3710  Sum=0. 

3720  ladd'O 

3730  IF  Ipos-I  THEN  Iadd-40 

3740! 

3750  !   CORRECTION  FOR  BROKEN  T/C  CHANNELS  5,  AND  47  CURVED  HEATER  ONLY. 

3760! 

3770  IF  Itype-1  THEN 

3780  FOR  1=0  TO  6 

3790  IF  (lpos=0  AND  1-5)  OR  (Ipos-1  AND  1-6)  THEN  3810 

3800  Sum=Sum+T( I+Iadd) 

3810  NEXT  I 

3820  Tijin>Su(n/6 

3830  ELSE 

3840  FOR  I-O  TO  7 

3850  IF  (Ipos-0  AND  1-7)  OR  (Ipos-t  AND  I-G)  THEN  3870 

3860  Sum-Sum+Tl I+Iadd) 

3870  NEXT  I 

3880  Tijin-Suin/7 

3890  END    IF 

3900  PRINT  USING  "12X.""Avg  Local  Heated  Wall  Temp  at  Inlet  (Tuiin)      -  "',3D. 

DD,""  (C) :Tuin 

3910! 

3920  !   CALCULATE  AVERAGE  LOCAL  HEATED  WALL  TEMPERATURE  NEAR  MIDDLE. 

3930! 

3940  Sum=0. 

3950! 

3950  !   CORRECTION  FOR  BROKEN  T/C  CHANNEL  15  CURVED  HEATER  ONLY,  AND  T/C 

3961!CHANNEL  11  STRAIGHT  HEATER  ONLY. 

3970! 

3980  IF  Itype-I  THEN 

3990  FOR  1-9  TO  15 

4000  IF  (lpos  =  0  AND  I-IO)  OR  (  lpos  =  0  AND  I- IS)  OR  (Ipos-1  AND  1  =  9)  THEN  4020 

401  0  Sum  =  Sum  +  Tl  I<-Iadd) 

4020  NEXT  I 

4030  IF  Ipos-0  THEN  Tumid-Sum/6 . 0 

4031  IF  Ipos-1  THEN  Tumid-Sum/ 7  .  0 
4040  ELSE 

4050  FOR  1-9  TO  16 

4060  IF  (lpos=0  AND  I-IO)  OR  (Ipos-1  AND  1=9)  OR  (Ipos-0  AND  I-ll)  THEN  4080 

4070  Sum=Sun+T( I+Iadd) 

4080  NEXT  I 

4081  IF  Ipos-0  THEN  Tuimid-Sum/6 . 0 
4090  IF  Ipos=1  THEN  Turn id=Sum/7 . 0 
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4100  END  IF 

4110  PRINT  USING  '•12X,""flvg  Local  Heated  Wall  Temp  Near  Middle  (Twmid)  -  "",30, 

DD.""  (C) ; Turn  id 

4120! 

4130  !   CALCULATE  A^^ERAGE  LOCAL  HEATED  UALL  TEMPERATURE  AT  OUTLET. 

4140! 

4150  Sum-0. 

4160! 

4170  !   CORRECTION  FOR  BROKEN  T/C  CHANNEL  26  CURVED  HEATER  ONLY. 

4180! 

4190  IF  Itype-1  AND  Ipos-0  THEN 

4200  FOR  1-23  TO  29 

4210  IF  I»26  THEN  4230 

4220  Sum  =  Sum  +  T(  I) 

4230  NEXT  I 

4240  Tuout=Sum/8 

42S0  ELSE 

4260  FOR  1-23  TO  29 

4270  Sum  =  Suni  +  T(I  +  Iadd) 

4280  NEXT  I 

4290  Tijout  =  Sum/7 

4300  END    IF 

4310  PRINT  USING  "12X,""Avg  Local  Heated  Wall  Temp  at  Outlet  (Twout)    -  ''■,3D. 

DD,""  (C) """  :Tu)out 

4320  PRINT 

4330  PRINT 

4340! 

4350  !   ENTER  ORIFICE  PRESSURE  DATA. 

4360! 

4370  BEEP 

4380  INPUT  "ENTER  PATM  (inHg),  DPM(inH20).  P1M(inH20).  RE  1"  .  Patm. Dpm,P1 m.Rel 

4390  PRINT  USING  "9X,""THE  FOLLOWING  ORIFICE  DATA  WERE  ENTERED:""" 

4400  PRINT 

4410  PRINT  USING  " 12X, " "Patml inHg)     DPM(inH2Q)     P1M(inH20)      RE1"" 

4420  PRINT  USING  " 1 4X,2(2D. DD, 1  OX ) , (2D.DD. 7X) ,5D. D" ;Patm.Dpm.P1m,Re1 

4430  PRINT 

4440! 

4450  !   CONVERT  PRESSURE  READINGS  TO  SI  UNITS  AND  P1M  TO  ABSOLUTE. 

4460! 

4470  Patm-=Patm*3376.8 

4480  Pdel=Dpm*248.84 

4490  P1=Patm-(P1m*248.84) 

4500  Torf=T(33) 

4510! 

4520  !   PHYSICAL  PROPERTIES  AND  CONSTANTS. 

4530! 

4540  R=286.987  !FOR  AIR. 

4550  Rho=Pl/(R*(Tort+273. 15) )  !DENSITY  OF  AIR  AT  ORIF  PLATE. 

4550  Gamma=1 .40  IFOR  AIR. 

4570  Cp=!006  ICp  GOOD  FOR  Tblk  BETWEEN  12  AND  33  (DEC  C). 

4580  Gc=1.0 

4590  Sigma=5.659E-8 

4600  Ecu=.12   !EMISSIVITY  OF  COPPER. 

4610! 

4620  !   DIMENSIONS  OF  CHANNEL. 

4630! 

4640  Dc-. 006350 

4650  Wid'.2540 

4650  Ri=.297 

4670  Pujet  =  2»(Dc-t-Wld) 

4680  Ac-Dc^Wid 
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4G90  Dhd  =  '^*Ac/Puet 

4700  Apl=. 0774192 

4710! 

4720  !   PROPERTIES  OF  INSULATION. 

4730! 

4740  Xins-.0127 

4750  Kins-3.8E-2 

4760! 

4770  !   ENTER  ORIFICE  CONFIGURATION. 

4780! 

4790  BEEP 

4800  INPUT  "SELECT  DIAMETER  OF  ORIFICE  (0-.5334,  1-1.0755  (inches) )" ,Si ze 

4810  IF  Size=0  THEN  Dorf-. 013548 

4820  IF  Size=t  THEN  Dorf-. 027318 

4830! 

4840  !   DIMENSIONS  OF  PIPE  AND  ORIFICE  PLATE. 

4850! 

4860  A-(PI*Dorf  2)/4 

4870  Dpipe-. 051772 

4880  Apipe-PI*Dpipe'2/4 

4890  Beta=Dorf /Dpipe 

4900  PRINT  USING  "  1 2X  .  "'Dorf  ( m )        A(ni"2)         Beta""'   . 

4910  PRINT  USING  " 1 2X , ( Z . 6D.6X ) , ( Z . 3DE ,6X ) , Z . 4D" ; Dorf , A . Be ta 

4920  PRINT 

4930  PRINT 

4940! 

4950  !   CORRELATION  FOR  EXPANSION  FACTOR  BASED  ON  ID  AND  t/2D  TAPS. 

4960  ! 

49  70  Y-l-(.41+.35*Beta'4)*(Pdel/(Gamma*Pl  )) 

4980! 

4990  !   CALCULATION  OF  FLOU  COEFFICIENT  (K)  IN  SI  UNITS. 

5000! 

5010  B-. 0002+2. 794E-5/DpipeM. 0038+1. 016E-5/Dpipe)*( Beta •2+(l 6. 5> I .968504E+2*Dp 

ipe)*Beta" 16) 

5020! 

5030  !   Kt  AND  K2  ARE  DUMMY  VARIABLES  TO  CALCULATE  Ko. 

5040! 

5050  Kl-.6014-5.3974064E-3/Dpipe-.25 

5060  K2=( .376+2.8971138E-2/Dpipe'.25)*(1.6129E-7/(Dpipe'2*Beta'2+6.35E-5*Dpipe) 

+Beta"4+1 .5*Beta" 16) 

5070  KO-K1+K2 

5080  K  =  Ko  +  100Q*B/ReT  .5 

5090! 

5100  I   Mu  CORRELATION  GOOD  FOR  Torf  BETWEEN  17  AND  41  DEGREES  CELSIUS. 

5110! 

5120  Mu=4.6971429E-8*Torf+I .7194  722E-5 

5130! 

5140  !   CALCULATE  Mdot  AND  REpipe  AND  COMPARE  TO  PREDICTED  REl. 

5150! 

5160  Mdot  =  K^A*r*(2*Gc*Rho*Pdei  )  '  .5 

5170  Repipe=(Mdot*Dpipe)/(Mu*Apipe) 

5180  Dif f-(ABS(Rel-Repipe)/Rel  )*iOO 

5190  Rel=Repipe 

5200  IF  Diff>.001  THEN  5080 

5210  PRINT  USING  "9X.""THE  FOLLOUING  DATA  WERE  CALCULATED:""" 

5220  PRINT 

5230  PRINT  USING  " 1 2X . " "Or i f ice  Expansion  Factor  (Y)     =   "",Z.4D,";Y 

5240  PRINT  USING  " 1 2X . " "Or  i  f  ice  Flow  Coefficient  (K)     =   "".Z.4D.":K 

5250  PRINT  USING  *  1 2X . " "Dens i ty  Based  on  Torf  ( Rho )      -   "",Z.4D,""     (ka/m' 

3) ;Rho 
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5260  PRINT  USING  " 12X, " "Viscosi ty  Based  on  Torf  (Mu)     -  "",2Z.3DE.""  (kg/m.s) 

;Mu 

5270! 

5280  !   Mu  CORRELATION  GOOD  FOR  Tblk  BETWEEN  T7  AND  41  DEGREES  CELSIUS. 

5290! 

5300  Mii-4.697t429E-8*Tblk  +  1  .7194722E-5 

5310  Red=(Mdot*Dc)/(Mu*Ac) 

5320  Rehd-(fidot*Dhd)/(Mu»Ac) 

5330! 

5340  !   CALCULATION  OF  POWER  INTO  HEATER  PLATE. 

5350! 

5360  Iadd-0 

5370  IF  Ipos>1  THEN  Iadd-39 

5380  IF  Ipos'O  THEN  Rpr-2. 01729 

5390  IF  Ipos-1  THEN  Rpr-2. 00839 

5400  Qp-(Emf (38+Iadd)*Emf (39+Iadd))/Rpr 

5410! 

5420  !   CALCULATION  OF  HEAT  INPUT  TO  AIR. 

5430! 

5440  Qair'Mdot*Cp*(Tout-Tin) 

5450! 

5460  !   CALCULATE  LOCAL  HEAT  TRANSFER  COEF.  AND  NUSSELT  NUMBER. 

5470! 

5480  Ka i r- 7. 7257 143E-5*Tcin+. 024 165836 

5490  Hin  =  Qair/(ApU(  Twin-Tcin)  ) 

5S00  Nuin-(Hin*Dhd)/Kair 

5510  Ka  i  r  =  7.  72571  43E-5*  Tern  id^-.  024 155835 

5520  Hmid-Qair/(Apl*(Tumid-Tcmid) ) 

5530  Numid=(Hmid»Dhd)/Kair 

5540  Ka  1  r- 7. 7257 143E-5^Tcout  +  . 024 165836 

5550  Hout=Qair/(Apl*(Tuout-Tcout) ) 

5560  Nuoijt-(Hout»Dhd)/Kair 

5570! 

5580  !   Kair  CORRELATION  GOOD  FOR  Tblk  BETWEEN  17  AND  41  DEGREES  CELSIUS. 

5590! 

5600  Ka 1 r- 7. 7257 143E-5* Tblk +.024 165336 

5610! 

5620  !   CALCULATE  AVERAGE  HEAT  TRANSFER  COEF.  AND  NUSSELT  NUMBER. 

5630! 

5640  Havg=Qair/(Apl*Tdel) 

5650  Nuavg-(Havg*Dhd)/Kair 

5660! 

5670  !   CALCULATE  HEAT  LOSSES. 

5580! 

5690  QloI-(Kins*Apl*(Tinsol-Tinso2) )/Xins 

5700  Qlo2-(Kins*Apl»( T inso2-Tinso3 ) )/Xins 

5710  QIo=(QIq!+Q1o2)/2 

5  720  QlilMKins*Apl#(  T  ins  1 1 -T  ins  i2  )  ) /X  i  ns 

5730  01  i2=(Kins-»rtpl  M  T  i  ns  i2- T  ins  i3  )  ) /X  i  ns 

5740  Qli=(Qlil+Qli2)/2 

5750  Rr  =  2.0*(C 1 -Ecu) / ( Ap MEcu ) )  +  ( I .0/Apl  ) 

5750  Qr  =  (Sigma*(ABS(  (  Tijo+273  .  1  5  )  '  4- (  Tu  i +273  .  1  5  )  "  4  )  )  ) /Rr 

5  7  70  Qdel-Qp-Qair-Qlol-Qr-Qli1 

5780! 

5790  !   CALCULATE  DEAN  NUMBER. 

5800! 

5810  De=.5*Red*( (Dc/2)/Ri)  '  .5 

5320  PRINT  USING  "  1 2X , '■  "Viscosi  ty  Based  on  Tblk  (Mu)     -  "".2Z.3DE,""  (kg/m.s) 

' ;  M  u 

5830  PRINT  USING  " 12X, " "Therm  Cond  Based  on  Tblk  (Kair)  -   "",Z.3DE,""  (W/m.K)" 
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" " :  K 

5840 
5850 
5860 
5870 
5880 
5890 
5900 
5910 
5920 
5930 
5940 
5950 
59bO 
5970 
atts 
5980 
(Uat 
5990 
(Uat 
6000 
(Wat 
6010 
(Wat 
6020 
6030 
in.N 
6040 
mid. 
6050 
ou  t , 
6060 
avg , 
6070 
6080 
6090 
6100 

61  10 
6120 
6130 
6140 
6150 
6160 
6170 
6180 
6190 
6200 
6210 
6220 

62  30 
6240 
6250 
6260 
6270 
6280 
6290 
6300 
6310 
6320 
6330 


a  1  r 
PRINT 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 


'12X.""P3tm 
'12X,""Pdel 

12X,""P1 
'12X,""Mdot 

1 2X , " "Rep ipe 
'12X.""Red 

12X.""Rehd 


.6D.DD." 
"".4D.DD 
,5D.DD." 
"".Z.4D 
,6D.D":R 
.6D.D";R 
.5D.D";R 


I 


ONLY  PRINT  DEAN  NUMBER  WHEN  HEATING 


IF  Itype-l  THEN  PRINT  USING  ■'12X.""De 

PRINT 

PRINT  USING  "12X.""Qp 
) ;  Qp  .  Qa  I  r 

PRINT  USING  '*12X.""Qlol 
ts) :Qlol .Oli1 

PRINT  USING  "12X,-"QIo2 
ts) ;Qlo2.Qli2 

PRINT  USING  ''12X.""Qlo 
ts)""";01o,Qli 


■Qr 


12X, 
12X, 


-Hin 
"Hmid 
"Hout 
"Havg 


.3D. 3D,-' 
■.2D. 3D.' 
",2D.3D.' 
•.2D. 3D,' 
'',2D.3D,' 


,3D.DD."' 
.3D.DD,'" 
,3D.DD."' 
,3D.DD."' 


PRINT  USING  "12X 
ts) :Qr.Qdel 

PRINT 

PRINT  USING  -12X 
ui  n 

PRINT  USING  "12X 
Num  id 

PRINT  USING 
Nuou  t 

PRINT  USING 
Nuavq 

PRiNT 

PRINT 

BEEP 

INPUT  "WILL  THERE  BE  ANOTHER  RUN?  (1-Y 

IF  Go_on<>0  THEN 

PRINT  USING  '"i,//- 

IF  lm=0  THEN  WAIT  300 

GOTO  1700 

END  IF 
'PRINT  USING  "8X.""END  OF  RUN""" 

PRINT  USING  "Li,*" 

PRINTER  IS  1 

ASSIGN  *File  TO  * 

END 

!   CONVERTS  EHF  TO  DEGREES  CELSIUS. 
j 

SUB  Tvsv(V.T) 

(:0M  /Co/  D(7)  .Aa(76.2) 

S  u  m  ^  0 

FOR  1=0  TO  7 

Sum=Sum+D( I )*V* I 

NEXT  I 

T  =  Sum 

SUBEND 
! 

!   CALIBRATES  T/C  READINGS. 


(N/m"2)""";Patm 
"  (N/m'2)""";Pdel 
(N/m-2)""";P1 
"  (kg/s)"""  ;ridot 

epipe 

ed 

ehd 

CURVED  SECTION. 

"".4D.D";De 
(Watts)"". 7X,"*Qair  - 
(Watts)"', 7X.""Qlil  • 
(Watts)"", 7X, ""0112  ■ 
(Watts)"", 7X.""Qli   ■ 
(Watts)"". 7X,""Qdel  ■ 

(W/m-2C)*".7X,""Nuin 
(W/m'2C)"",7X,""Numid 
(W/m-2C)"",7X.""Nuout 
(W/m-2C)"",7X,""Nuavg 


,0-N)",Go_on 


,3D.3D,""  (H 

"",2D.3D,"" 

•,2D.3D,"" 

".2D.3D,'" 

-.2D. 3D,"" 

".3D.DD.";H 
",3D.DD,";H 
".3D.DD,";H 
",3D.DD,";H 
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63^0! 

6350  DEF  FNTcorr(T.I) 

6360  COM  /Co/  D(7) ,Aa(76,2) 

6370  Tc=Aa(I.O) 

6380  FOR  J=1  TO  2 

6390  Tc-rc>Aa(I.J)*T*J 

6400  NEXT  J 

6410  RETURN  Tc 

6420  FNEND 
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APPENDIX  B 
SAMPLE  CALCULATION 

Figure  (B.l)  shows  the  components  of  the  heat  balance  in  the  test  channel.  A 
sample  calculation  follows  that  details  the  methods  used  in  the  data  reduction  program 
to  estimate  the  value  of  the  energy  balance  components  and  then  calculate  the 
hydraulic  Reynolds  number  and  average  Nusselt  number  for  each  data  set.  A  copy  of 
an  actual  record  of  one  set  of  experimental  data  for  the  concave  test  section  precedes 
the  sample  calculations  together  with  a  list  of  the  parameters  necessar>'  to  make  the 
sample  calculation. 


Qp    — 

r—HEATED    COPPER    WALL 

..J\^ 

i                                   '     [ 

'    1°'             "-^0.:, 

\     1 

-? 

X 

.._..._;_._.. 

Figure  B.l     Test  Section  Energy  Balance. 
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1.       SAMPLE  CALCULATION  DATA 


Month,  Day,  and  Time:  09 : 1 b : 1 5 : 1 4 : 54 

THE  FOLLOWING  DATA  WERE  RECORDED: 

Heater  Uoltaae  (Vh)       =  91.7E.B  (V) 
Resistor  Voltage  i\iT)  =   7.704  (V) 

Precision  Resistor  (Rpri  =   1'.  01  729  lOhms ) 

OUTER  PLATE: 

T/C  flumoer:      12  3       4       5       b 

Temp  lO:      72.33    63.0  7  7b.  00  &8.01    7&.30    72.90 

T/C  Nufnt.er :      7       8  9      10       I  1       12 

Temp  (C):     79.99   72.13  73.77  74.26   77. U3   7B.00 

T/C  Number:     13      14  IE.      16      17       18 

Temp  (O:     74.32   7B.52  7B.70  53.57   62.45   78.73 

T/C  Number:     19      20  21      22      23      24 

Temp  (C):     79.13   84.26  85.87  69.58   78.27    78.07 

T/C  Number:     25      26  2  7      28      29      30 

Temp  (O:      77.53   80.33  47.70  79.25    80.06    83.47 

OUTER  INSULATION: 

T/C  Number:     31  (Tinsol)  32  (Tinso2)    33  (Tinso3) 

Temp  (O:         82.38  55.88         40.01 

ORIFICE  TEMP  (Torf)  -   30.25  (C) 

INLET  TEMPERATURE: 

T/C  Number:     35      35  37      38 

Temp  (C):     22.09   22.07  22.08  22.07 

INNER  PLATE: 

T/C  Number:     41       42  43       44       45       45 

Temp  (Cj:     23.89   23.69  23.40  23.34   23.21    23.49 

T/C  Number:     47      48  49      50      51       52 

Temp  (O:      23.28   26.29  23.55  23.27   24.03    24.23 

T/C  Number:    53      54  55      56      57      58 

Temp  (O:     23.71    23.82  23.59  23.53   23.56   24.02 

T/C  Number:     59      60  61      62      63      64 

Temp  (C):      23.80    24.49  24.27  24.87    24.68    25.87 

T/C  Number:     65      6b  67      68      59      70 

Temp  (O:     25.56   25.32  25.15  25.04   24.95   24.93 

INNER  INSULATION: 

T/C  Number:     71  (Tinsil)  72  (Tinsi2)    73  (Tinsi3) 

Temp  (O:         24.40  24.13         23.70 

OUTLET  TEMPERATURE: 

T/C  Number:     74      75  75      7? 

Temp  (O:      30.55    30.73  30.65  30.59 
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THE  FOLLOWING  TEMPERATURES  WERE  CALCULATED: 

Average  Outer  Wall  Temperature  (Two) 

Average  Inner  Wail  Temperature  (Tui) 

Average  Outlet  Temperature  (Tout) 

Average  Inlet  Temperature  (Tin) 

Channel  Inlet  and  Outlet  Temp  Diff  (Tdiff) 

Average  Bulk  Temperature  (TbU) 

Mean  Temperature  Difference  (Tdel) 

Outer  and  Inner  Wall  Temp  Difference  (Tudiff) 

Average  Local  Channel  Temp  at  Inlet  (Tcin) 

Average  Local  Channel  Temp  at  Middle  (Tcmid* 

Average  Local  Channel  Temp  at  Outlet  (Tcout) 

Avg  Local  Heated  Wall  Temp  at  Inlet  (Twin) 

Avg  Local  Heated  Wall  Temp  Near  Middle  (Twmidl 

Avg  Local  Heated  Wall  Temp  at  Outlet  (Tuout) 


THE  FOLLOWING  ORIFICE  DATA  WERE  ENTERED: 

Patm(inHg)      DPM(inH20)  PlM(inH20)       REl 
29.86          17.60  16.90       15000.0 

Dorf(m)        A(m'2)         Beta 

0.027318      5.361E-0A  0.S277 


THE  FOLLOWING  DATA  WERE  CALCULATED: 

Orifice  Expansion  Factor  (Y) 
Orifice  Flow  Coefficient  (K) 
Density  Based  on  Torf  (Rho) 
Viscosity  Based  on  Torf  (Mu) 
Viscosity  Based  on  TbU  (Mu) 
Therm  Cond  Based  on  Tblk  (Kair) 

Patm    =■    100848.13  (N/m"2) 


76, 

.46 

(C) 

24. 

,  16 

(C) 

30, 

,68 

(C) 

-lO 

,08 

(C) 

^s'. 

.61 

(C) 

26. 

.3S 

(C) 

50, 

,08 

(C) 

52. 

,30 

(C) 

oo 

,44 

(C) 

25! 

,30 

(C) 

30, 

,32 

(C) 

73. 

.62 

(C) 

73, 

,06 

(C) 

79. 

,78 

(C) 

0 

.9859 

0. 

,6326 

1 , 

.1099 

(kg/m"3) 

18, 

,6I6E- 

■06 

(kg/m.s) 

18 

.434E- 

-06 

(  k  q  /  m  .  3  ) 

2, 

,b20E- 

■02 

(W/m.K) 

Pdel 

4379.58  (N/m-2) 

PI 

=   96642.74  (N/m'2) 

lidot 

0.0360  (ka/s) 

Rep  ipe 

=   47611.8 

Red 

7697.2 

Rehd 

=   15018.9 

De 

397.9 

Qp 

> 

350.411  (Watts) 

i:  1 0 1 

= 

5.908  (Watts) 

0io2 

= 

3.908  (Watts) 

Qlo 

= 

4.303  (Watts) 

Or 

= 

1.396  (Watts) 

Hin 

= 

73.7^  (W/m'2C) 

Mm  id 

= 

34.33  iW/m-2C) 

Hout 

= 

81  .48  (W/m"2C) 

Havg 

= 

30.^7  (W/m-2C) 

Qair  = 

31  I .991 

0  1  1  1  = 

.063 

131  i2  = 

.  100 

Oil   = 

.031 

IJdel  = 

30.453 

flu  in 

37.67 

Humid  = 

40.03 

Nuout  = 

38.08 

Nuavg  = 

38.05 

(Watts) 

(Watts) 
(Watts) 
I Wa 1 1  3  ) 
( Uatts) 


END  OF  RUN 
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V,  =  91.755  V  P^.^^  =. 5207  m 

Vp^  =  7.704  V  rj  =  .297  m 

Rp^  =  2.01729  n 

Patm  =  29.86  in.  Hg 

AP  =   17.60  in.  H^O 

Pj  =  16.90  in.  H^O 

Tj^  =  22.08  °C 

Tout  =  30.68  X 

\o  =  ^6.46  X 

\,,  =24.I6°C 

Tinsil  =  24.40  °C 

Tins.2  =24.13"C 

Tinsi3  =  23.70  =C 

Tinsel  =  82.38  X 

Tinso2  =  56.88  °C 

Tinso3  =  ^0.01  X 

Torf  =  30.25  °C 

K^j^  =  .0262  W,  m-°K 

Cp  =   1007JKg-°K 

^*air  =  ^-6971  X  10-8  X  T^.^  +  1.71947  X  10-5    j.g,j^.5 

^ins  =  -0380  W/m-X 

AX.„3  =.01 27  m 

^cu  =  012 

ff  =  5.669  X  10  8  W;m^-°K'* 

P  =  .5277 

Y  =  1.40 


g-  =1.0  Kg-m/N-s 


2 


-c 

R  =  286.98  Xm,  Kg-°K 

F  =1.0 

wo-vvi 

Ap,  =  .07742  m^ 

A^  =  .0016  nr 

A  .  =  .00211  nr 

pipe  -  ^^i  I  111 

D^  =  .00635  m 

D  .  =  .0518  m 
pipe  .V-  ^  lu  in 

D^^^  =  .027318  m 
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2.  TEMPERATURE  CALCULATIONS 

a.  Average  Temperature  of  the  Outer  (Inner)  Wall 

V  Outer  (Inner)  Plate  Temperatures 

T     (T  ..)  =  —^ ^ =  76.46''C  (24.16°C) 

wo     wi        Number  of  Outer  (Inner)  Plate  Thermocouples 

b.  Average  Channel  Outlet  (Inlet)  Temperature 

V  Outlet  (Inlet)  Temperatures 

T    ,(T  )  =  —=^ =  30.68X  (22.08°C) 

out     in        Number  of  Outlet  (Inlet)  Thermocouples 

c.  Channel  Inlet  to  Outlet  Temperature  Difference 

T.o-  =  T      -  T     =  30.68X  -  22.08X  =  8.6°C 

dill  out         in 

d.  Bulk.  Temperature  of  the  Flow  Through  the  Channel 

T.     +  T    ,  22.0S^C  +  30.68X 

T. ,.    =  -^ ^i^    =    =    26.38°C 

^"^  2.0  2.0 

e.  Average  Temperature  Difference  of  the  Heated  Wall  and  Fluid  Bulk  Temperature 

AT  =  T      -  T. ,.    =  76.46X  -  26.3SX  =  50.08X 
wo        blk 

f.  Temperature  Difference  Between  the  Outer  and  Inner  Walls 

T    ..^  =  T      -  T  .  =  76.46'^C  -  24.16X  =  52.30^C 
vvdin  wo         vvi 

3.  POWER  CALCULATIONS 

a.  Power  supplied  by  the  Heated  Plate 

V,  XV  91.755V  X  7.704V 

Q    =  -^ ^  =  =  350.41  Watts 

P  R  2.01729n 

pr 
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b.  Poner  Loss  Through  the  Outer  Wall  Insulation 

Qio  =  (i''2XQ,oi  +  Q102) 


T         -  T 

insol         inso2 


=  (i/2)[(-ins^i inio^)k.    A  , 

AX  '"^    P 

^  ins 


T         -  T 

_l_   /    inso2         inso3    vi.       a     1   = 
^  AX  ^  •ns-'^plJ 

ins 

.  (0.038  W  mX) 

(l/2)(.07742m^) [(82.38  -  56.88(°C))  +  (56.88  -  40.01(°C))]  =  4.908W 

.0127m 

c.  Power  Loss  Through  the  Inner  Wall  Insulation 

Q„  =  (1.2KQ„,  +  Q„,) 

=  (l;2)l(i^^-im2.)k.„,Ap, 

ins 

T        -  T 
^         AX.  '"'    P^ 

ins 

.  (0.038  W-mX) 

(l/2)(.07742m-)  [(24.40  -  24.13(^C))  +  24.13  -  23.70(^C))]  =  .08 IW 

.0127m 

d.  Heat  Loss  by  Radiation 

/.    Radiation  Resistance 

1  -  £  1 

R    =  2x( ^i^)  + 


A  ,  X  £  A     X  F 

pi         cu  pi          wo  VVl 

12  12                                 , 

( 1)  =  — -,(- 1)  =  202.36  m  ^ 


Qr  = 


A  ,        £  0.07742m-    0.12 

pi  cu 

2.    Heat  Radiated 


(5.669  X  10  ^)[(76.46+273)^  -  (24.16+273)^1 

^ ■ —  =   1.996  W 

202.36 
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4.        MASS  FLOW  RATE  CALCULATIONS 

a.  Pressure  Conversions 

P  ,     =  29.865  m  Hg  x  3376.8  =   100848.13  Nm^ 
aim  ° 

\?  =  17.60  in  H.O  x  248.84  =  4379.58  N  nr 

P,  =  (16.90  in  H^O  X  248.84)-  100848.13  =  96642.74  N;m^ 


b.  Density  of  Air 


P,                        (96642.74)  , 

p  .,  =  — i —  =  =   1.1099  Kg,m^ 


RT^^f.        (286.98)  (30.25  +  273) 
c.   Expansion  Factor 
Y  =   1  -  (0.41  +  0.35p-*) 


VPi 


4379.58 
=   1  -  [0.41  +  0.35(0.5277)-^] — =  0.9859 


(1.40)(96642.74) 


d.  Mass  Flow  Rate 


m  =  YKA  v'2gcPair^P  =  (.9859)K(.0005861m^) 

X  V2(l  Kg-m,Ns-)(  1.1099  Kg  nr )(4379.58  N/m^)  =  .05669 K  Kg/s 

Iterate  to  find  m: 

Assume  a  hydraulic  Reynolds  number  =  15000 

D  .    A  .0518  X. 0016 

Re.      -  (     P'P^    ^   )  Re.  ,  =  ( )-  3.1961  x  Re.  , 

P'P«  D,.A./      hd       \01229x.002ir  ^^^ 

hd     pipe 

-♦  Re  .      -  15000(3.1961)  =  47942 
pipe  ^ 

Interpolate  in  the  appropriate  table  in  Reference  15  for  K,  the  How  coetricient: 
K  =  .6329 

{Note   [hat   [he    computer  program   used   to    reduce    the   data   used  an   approximating 
polynomial  to  find  the  correct  value  of  the  Jlow  coejficient,  K.) 

->  m  =  (.05669)(.6329)  =  .0359  Kg  s 

Solve  for  the  improved  value  of  Re  .    : 
^  pipe 

m  D  .                    (.0359K.O518) 
Re  .      =  2121-  =  ^ ^ =  47315 

^'^'      ''"Vpe^^air  (.0021  1 )( 18.616  X  10  ^) 

Iterate  until  successive  values  of  m  dilfer  by  less  than  .001:  rn  =  .0360 
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REYNOLDS  NUMBER  CALCULATIONS 


a.  Re  Based  on  Pipe  Diameter 

niD.  (.0360)(.0518) 

Re  .      =  EiE^  =  ^ -^  =  47475 

P'P^      A  .    u  .         (.00211X18.616  X  10^) 

b.  Re  Based  on  Channel  Height 

m  D  (.0360)(.00635) 

Rg     =:  C        _  ^ ^ '  =   7751 

'^        A^H^jj.  (.0016X18.434X10"^) 

c.  Re  Based  on  Hydraulic  Diameter  of  the  Channel 

rhD,.  (.0360X.01229) 

Re.  .  =  ^^  =  ^ ^ -—  =   15000.8 

^^         Au.  (.0016H18.434X  10  ^) 


6.        HEAT  TRANSFERRED  TO  THE  WORKING  FLUID 

Q^;^  =  mcp  (T^^^  -  T.J  =  (.0360  Kg,  sX  1007  J/Kg-°CX30.86  -  22.08rC  =  312  W 


7.        AVERAGE  HEAT  TRANSFER  COEFFICIENT  CALCULATION 

_        Q.  312  W  , 

h  =  -laii_  = =  80.5  W;m^-°C 

A  jAT  (.07742m^X50.08°C) 


8.        AVERAGE  NUSSELT  NUMBER  CALCULATION 

_  h  D.  .         (80.5  W/m'-°CX.01229m-) 

^Ah  =  ^^^—  = =  37.8 

^^        K  .  .0262  W.m-'C 

air 


9.        DEAN  NUMBER  CALCULATION 


Re.  /D  ;2  7751         /. 00635/ 2 

De  =  — ^  =  V -^ —  =  y =  400.7 

2  r.  2  .297 

1 
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APPENDIX  C 
EXPERIMENTAL  UNCERTAINTY 

The  experimental  uncertainty  in  calculated  parameters  was  determined  in 
accordance  with  the  method  described  by  Kline  and  McClintock  [Ref.  23]. 

The  qualitative  results  shown  here  correspond  to  the  data  set  and  calculation  of 
Appendix  B.  There  was  a  significant  decrease  in  the  experimental  uncertainty  of  the 
average  Nusselt  number  for  Reynolds  numbers  above  17,000  in  this  experimentation 
compared  to  that  of  Galyo  [Ref  14:  pp.  104-105]  because  of  the  changes  made  to  the 
test  channel  that  promoted  mixing  of  the  working  fluid  at  the  exit  of  the  channel. 

The  equations  used  to  determine  the  experimental  uncertainty  of  calculated 
parameters  are  shown  below. 

1.  DENSITY  OF  THE  WORKING  FLUID 

8p  6R     ^  8P,    -,  8T    f     ->  n  ? 

—  =  [(  -r-  )'    +    {  -^  )'    +    (  -r"^  9  f-' 

P  R  Pi  T^rf 

2.  MASS  FLOW  RATE  OF  THE  WORKING  FLUID 

6m  8Y     -,  6K    T  6 A   -, 

^  =  [((  -r-  )-    +    (  -^  )^    +    (  — -)-  } 
m  Y  k  A 


6p  .         ,  SAP    ,      r,  - 


p  .  A? 

•  air 


3.        HEAT  TRANSFERRED  TO  THE  WORKING  FLUID 

§0  .  6m     ^  6c       ^ 

-7^  =  [{  -r-  )-   +   (  —^  r 


(T     ,  -  T.    ) 
^    out         in 


)- 
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4.   AVERAGE  HEAT  TRANSFER  COEFFICIENT 

6h  6Q,.       .  8A  .    ,  SAT  ,     .  , 

h  Qair  Apl  ^T 


5.        AVERAGE  NUSSELT  NUMBER 


Tid     _   r/  n2     ^     /  hd  \2     _^     /  air.  \2    i  0.5 


>^A.d  ^  c>hd  K 


air 


6.        HYDRAULIC  REYNOLDS  NUMBER 

8Re  6m     ,  8D.  .  .,  6ii  .        ,  5A      ->     ^  ,. 

Re  ^  ^hd  ^*air  ^c 


Several  parameters  used  in  the  calculations  described  herein  had  uncertainties 
that  were  constant  over  the  entire  range  of  testing.  These  quantities  and  their  relative 
uncertainties  follow. 

Quantity  (B)  Uncertainty  (  ±  6B  B) 


A  .0020 

A  .0408 

c 

A  ,  .0006 

pl 

c  .0020 

P 

D^  .0400 

Dhd  -0391 

K  .0030 

K  .  .0038 

air 

aim 

R  .0003 

W. .  .0080 

id 

Y  .0020 

u  .  .0054 

•^air 
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The  approximating  polynomials  used  to  resolve  thermocouple  voltages  into 
temperatures  matched  actual  temperatures  to  within  less  thaniO.rC  over  a  range  of 
temperatures  from  20°C  to  85X;  however,  in  some  cases  this  temperature  range  was 
exceeded.  A  conservative  estimate  of  the  uncertainty  of  temperatures  measured  by 
thermocouples  was  set  at  ±0.rC,  the  accuracy  of  the  thermocouples  specified  by  the 
manufacturer  of  the  thermocouple  wire.  The  following  uncertainties  were  used  in 
calculations: 

Quantity  (T)  Uncertainty  (  ±  T) 


T 

T. 

in 

T 

orf 

out 

T     u 

\vh 

T       -  T. 

out         in 

AT 


0.1 

X 

O.I 

°c 

O.I 

°c 

O.I 

°c 

0.1 

X 

0.2 

°c 

0.2 

°c 

For  the  data  taken  at  a  hydraulic  Reynolds  number  of  15000  in  the  concave  test 
section  the  following  additional  uncertainties  applied: 

Quantity  (P)  Uncertainty  (  ±  P) 


Pj  0.05 

AP  O.IO 


The  results  of  the  uncertainty  calculations  at  a  Reynolds  number  of  15000  in  the 
concave  test  section  were: 

Quantity  (B)  Uncertainty  (  ±6B;B) 


p  .00445 

rn  .00544 

Qair  -02397 

h  .02431 

Nu^j  .04620 

Rej^j  .05703 


The  major  contribution  to  experimental  uncertainty  in  both  the  average  Nusselt 
number  and  hydraulic  Reynolds  number  was  the  uncertainty  in  cross  sectional  area  of 
the  channel.  The  hydraulic  diameter  of  the  channel  was  a  function  of  channel  cross 
sectional  area,  resulting  in  a  relatively  large  uncertainty  in  the  hydraulic  diameter. 
These  two  parameters  accounted  for  approximately  fourty  percent  of  the  uncertainty  in 
the  average  Nusselt  number  calculation,  and  about  eighty  percent  of  the  uncertainty  in 
the  calculated  hvdraulic  Revnolds  number. 
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